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ABSTRACT 

An  analytic  model  is  developed  to  allow  derivation  of  the 
probability  that  a  molecule  in  the  atmosphere  of  Saturn's 
rings  collides  with  at  least  one  ring  particle  when  tra- 
versing the  ring  plane.   The  resulting  expression  involves 
details  of  the  molecule's  trajectory,  including  the  velocity 
of  the  molecule  relative  to  ring  particles;  thus,  the  theory 
is  compatible  with  recently  developed  ballistic-transport 
computer  models  used  in  the  analysis  of  the  rings.   The 
collision  theory  is  applied  to  the  case  of  a  low  energy, 
isotropic  molecular  production  process  to  make  inferences  on 
the  spatial  extent  of  the  ring  atmosphere  associated  with 
such  a  source.   The  high  frequency  of  collisions  expected 
for  the  atmosphere  in  the  vicinity  of  the  A  and  B  rings 
suggests  a  toroidal  atmosphere. 
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LIST  OF  SYMBOLS 


b         thickness  of  a  sublayer  of  a  "slab"  whose  overall 
thickness  is  h  (h  <  b) 

h         thickness  of  the  translating  slab  (containing  par- 
ticles) through  which  a  molecule  passes  (corres- 
ponds to  ring  thickness  in  application  of  the 
collision  theory) 

n(r)       number  density  of  particles  in  a  slab  (or  in  the 
ring  for  applications  of  the  collision  theory) 

r         particle  radius  (sometimes  referred  to  in  this 
paper  as  the  "size"  of  a  particle) 

v         ejection  velocity  of  a  molecule,  relative  to  its 
-*        source  particle 

v         velocity  of  a  molecule 

m  J 

-> 

v         velocity  with  which  the  slab  translates  (equivalent 

to  a  representative)  velocity  of  the  particles 

"contained"  in  the  slab 

D         distance  that  the  slab  translates  in  time  t 

L         scale  length  for  the  extent  of  the  slab  in  direc- 
tions normal  to  its  thickness 

N  equilibrium  number  of  molecules  (of  a  given  type) 

in  ring  atmosphere  deriving  from  a  specific  produc- 
tion mechanism 

N(r.)      number  of  slab  (ring)  particles  of  size  r-j_,  con- 
tained in  the  slab  (meaningful  only  for  discretely 
distributed  r) 

P  probability  that  a  molecule  collides  with  at  least 

one  ring  particle  during  a  given  passage  through 
the  slab  (ring) 

P   (r.)    probability  that  a  molecule  passes  through  the 

slab  (ring)  without  colliding  with  a  particle  of 
size  r.  (only  strictly  meaningful  for  discrete  r) 

P   ,.v     nrobability  that  j  particles  of  size  r-  are  not  in 
no  h      *_,      i  i  •  •        •  i 

J      the  collision  region 


Q  (r,A)    extinction  cross  section,  sometimes  referred  to  as 
the  extinction  coefficient  (relates  to  the  effec- 
tive size  of  a  particle  with  regard  to  removal  of 
electromagnetic  energy  from  a  beam  via  scattering 
and  absorption) 

R         orbital  radius  (see  Section  III) 

R         interaction  radius  (see  Appendix  B) 

S         sticking  coefficent  (the  probability  that  a  mole- 
cule will  adhere  to  a  surface  when  it  collides 
with  that  surface) 

T         a  "mean"  orbital  period  typical  of  the  trajectories 
of  all  molecules  of  a  certain  type  produced  by  a 
given  process 

T         number  of  completed  orbital  periods 

W  production  rate  (of  molecules)  associated  with  a 
specific  mechanism,  usually  including  only  those 
molecules  which  assume  closed  orbits 

a,6,Y     direction  cosines  which  specify  a  molecule's  path 
as  it  passes  through  the  slab,  when  working  in 
a  frame  of  reference  in  which  the  slab  translates 

at  v   (see  Figure  2) 

P        * 

a'  ,8' ,y'   direction  cosines  which  specify  a  molecule's  path 

through  the  slab,  when  working  in  a  frame  of  refer- 
ence in  which  the  slab  is  fixed  (see  Figure  C2) 

<t>  angle  used  to  describe  the  deviation  of  the  actual 

direction  of  v   from  the  assumption  that  v   is 

P  P 

constant      e  * 

X  wavelength  (of  electromagnetic  radiation) 

y         angle  that  an  electromagnetic  beam  makes  with  the 
surface  normal  of  a  slab  (ring) 

9         angle  in  the  x-y  plane  between  the  projection  of 
the  molecular  path  and  the  +x  direction  (vp 
direction) ,  defined  in  a  frame  of  reference  in 

-y 

which  v_  ^  0 ,  referred  to  in  this  paper  as  the 
"orbital  aximuth" 

9'         same  as  9,  except  defined  in  the  frame  of  reference 

with  v   =  0 
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I.   INTRODUCTION 

It  has  been  well  established  that  the  particles  of  which 
Saturn's  rings  are  composed  are  primarily  water  ice  (Cuzzi 
et  al .  ,  1984).   Accepting  this,  one  would  expect  that  various 
erosive  processes  would  produce  a  ring  atmosphere  comprised 
of  H20  molecules  and  their  constituents,  e.g.,  H  atoms.   The 
density  of  such  an  atmosphere  would  be  determined  by  balancing 
the  production  rates  (of  atoms  and  molecules)  associated  with 
the  erosion  mechanisms  with  losses  resulting  from  capture  by 
Saturn,  escape  from  the  Saturnian  system,  and  recapture  by  the 
ring  particles. 

Early  theoretical  estimates  of  a  ring  atmosphere  (Denne- 
feld,  1974;  Blamont,  1974)  considered  sources  such  as  subli- 
mation, meteroid  impact,  and  bombardment  by  solar  and 
interstellar  wind.    These  analyses  were  concerned  only  with 
the  atmosphere  associated  with  the  A  and  B  rings  (see  Figure  1) , 
by  far  the  most  significant  potential  sources  of  material  for 
a  ring  atmosphere,  owing  to  their  large  masses  (Appendix  A 
gives  a  summary  of  ring  nomenclature  and  dimensions) .   Denne- 

feld  calculated  the  density  of  the  H„0  atmosphere  to  be  "25 

-3 
cm   ,  the  molecules  confined  to  a  toroidal  region  closely 

surrounding  the  A  and  B  rings.   He  estimated  the  H  atom 


Dennefeld  presented  a  conceptual  model  for  the  study  of 
ring  atmospheres,  which  is  of  fundamental  importance.   This 
model  is  discussed  later  in  this  section. 


B       A     F 

Cassini  Division 


extends  to 


Figure  1.   Nomenclature  of  the  Saturnian  Ring  System 
(drawn  to  scale) 


component  of  the  atmosphere  to  have  a  density  of  -1  cm 

5 
distributed  through  a  sphere  of  ~3  x 10   km,  well  beyond  the 

outside  edge  of  the  A  ring. 

Schemes  for  detection  of  the  atmosphere  all  rely  on 

measurement  of  H  Lyman-a  emission  from  the  vicinity  of  the 

rings,  the  most  sensitive  means  available.   Thus,  analysis  is 

limited  to  the  neutral  H  atmosphere.   Using  this  approach, 

several  researchers  put  Dennef eld's  estimates  to  the  test. 

Observations  from  instruments  onboard  rockets  (Weiser  et  al . , 

1977),  satellites  (Barker,  et  al . ,  1980;  Clarke,  1981), 
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Pioneer  II  (Judge  et  al.,  1980) ,  and  Voyager  I  (Broadfoot  et 
al.,  1981)  detected  H  atoms  near  the  A  and  B  rings.   Results 
were  reasonably  consistent,  suggesting  a  H  density  2-3  orders 
of  magnitude  greater  than  Dennefeld  predicted.   The  shape 
of  the  atmosphere  (spherical  versus  toroidal)  cannot  be  con- 
clusively deduced  from  these  Lyman-a  measurements,  but  pre- 
liminary analysis  of  Voyager  data  indicates  that  it  does  not 
extend  beyond  the  outer  edge  of  the  A  ring  (Cuzzi  et  al .  , 
1984)  . 

Efforts  to  reconcile  theoretical  estimates  of  the  H 
atmosphere  with  observations  have  primarily  been  directed 
toward  proposing  new  sources  or  re-evaluating  those  previously 
considered.   Proton  flux  from  Saturn's  ionosphere  (Ip,  1978), 
photodissociation  of  H^O  on  the  ring  particles  (Carlson,  1980), 
magnetospheric  ion  sputtering  (Cheng  et  al . ,  1978,1980),  and 

meteroid  impact  (Morfill  et  al . ,  1983;  Ip,  1983)  are  among 

2 
the  most  seriously  studied  erosive  mechanisms.    Of  these 

meteroid  impact  appears  likely  to  be  the  dominant  source. 

While  understanding  of  the  ring  atmosphere  has  increased, 

Dennefeld' s  conceptual  model  has  remained  the  foundation  for 

its  theoretical  analysis.   It  is  worthwhile  to  review  this 

model.   Suppose  that  a  particular  erosive  process  produces 

atmospheric  molecules  (e.g.,  H„0  or  H)  with  these  molecules 


2 

Some  have  also  considered  H  atom  sources  external  to  the 

rings  such  as  atoms  escaping  from  Saturn's  atmosphere 
(Shemansky  and  Smith,  1982).   Preliminary  calculations  indi- 
cate that  this  process  may  be  quite  significant. 
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having  some  initial  velocity  distribution.   The  trajectory 

of  each  molecule  is  primarily  governed  by  its  gravitational 

interaction  with  Saturn,  the  dominant  central  body  in  the 

Saturnian  system.   If  the  ring  atmosphere  is  quite  tenuous, 

collisions  between  molecules  as  they  move  in  their  Keplerian 

orbits  is  unlikely  (this  assumption  poses  no  problem  for 

Saturn's  rings,  since  for  a  density  of  even  1000  cm    the 

mean  free  path  is  -10   km).   A  molecule's  initial  velocity 

and  position  within  the  ring  determines  its  orbit  for  which 

there  are  three  possible  cases  : 

(1)   the  trajectory  takes  the  molecule  into  Saturn's 

atmosphere,  thus  removing  it  from  the  ring  atmosphere 

(2).   the  orbit  is  open  (hyperbolic)  with  the  molecule 
escaping  the  entire  Saturnian  system. 

(3)   the  orbit  is  elliptic  with  the  molecule  passing 

through  the  ring  once  or  twice  during  each  orbital 
period  (depending  on  the  eccentricity  of  the  orbit) . 

For  those  molecules  with  elliptic  orbits,  some  will  be  lost 

from  the  atmosphere  to  the  rings.   The  probability  that  this 

occurs  during  each  passage  of  a  molecule  through  the  rings 

is  the  product  of  the  probabilities  of  its  colliding  with 

ring  particles  and  of  its  sticking  to  a  particle  during  this 

process.   The  spatial  extent  of  the  atmosphere  tends  to  be 

toroidal  if  the  initial  speed  of  the  molecules  relative  to 

the  orbital  speed  of  the  ring  particles  is  small.   On  the 


The  first  two  cases  apply  if  the  erosive  mechanism 
produces  molecules  with  high  speeds  relative  to  the  ring 
particles;  the  third  is  dominant  for  relatively  lower  energy 
processes . 
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other  hand  if  the  initial  speed  is  much  greater  than  that  of 
the  ring  particles,  then  the  atmosphere  will  tend  to  be  more 
spherical.  A  mathematical  expression  of  this  process  of  the 
form  introduced  by  Dennefeld  is 


N   = 


WT 

SP   ' 
c 


where 


N   =   number  of  molecules  (of  a  given  "kind")  in 
the  equilibrium  atmosphere  (from  a  specific 
production  mechanism) . 

W   =   production  rate  (sec   )  of  molecules  asso- 
ciated with  a  specific  mechanism  which 
assume   elliptic  orbits  (call  this  the 
"source  term" ) . 

T  e      mean  orbital  period  (sec);  this  term  incor- 
porates initial  velocity  information,  so  it 
might  more  generally  be  referred  to  as  the 
"ballistic  term." 

S   =   sticking  coefficient,  the  probability  that 
a  molecule  will  adhere  to  a  ring  particle 
when  it  collides  with  it. 

P    =   probability  that  the  molecule  collides  with 
a  ring  particle  as  it  passes  through  the 
ring. 3 

This  model  is  admittedly  crude,  but  it  does  contain  the 

essence  of  the  physical  system.   To  see  how  the  theory  has 

evolved,  each  term  is  discussed  below. 


4 

A  tacit  assumption  in  this  formulation  is  that  the 

average  trajectory  intersects  the  ring  only  once  during 
an  orbital  period.   If  it  passes  through  the  ring  twice 
during  a  period,  Pc  has  a  slightly  different  interpretation, 
but  the  other  terms  are  unaffected. 
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The  source  term,  W,  has  received  much  attention  as  dis- 
cussed in  some  detail  earlier. 

The  ballistic  term,  T,  is  receiving  increased  emphasis. 
Ballistic  transport  models  (Durisen,  1984;  Ip,  1983)  incor- 
porate inhomogeneities  found  in  the  macrostructure  of  the 
rings  with  improved  estimates  of  initial  velocity  distribu- 
tions to  allow  better  specification  of  molecular  trajectories. 

The  sticking  coefficient,  S,  is  probably  the  least  under- 
stood factor  in  the  theory.   The  physical  chemistry  under- 
lying the  phenomenon  is  not  well  developed,  while  laboratory 
experiments- do  not  appear  to  accurately  replicate  conditions 
found,  in  the  ring  environment.    While  it  may  be  possible  that 
the  physical  chemistry  problem  will  be  solved,  a  more  basic 
dilemma  will  likely  remain  for  quite  some  time:   S  cannot  be 
determined  accurately  unless  detailed  knowledge  of  individual 
ring  particles  is  developed.   Specifically,  impurity  levels 
and  the  phase  of  the  water  ice  in  the  particles  are  poorly 
known  (Weidenschilling ,  19  84),  and  the  surface  texture  of 
the  particles  is  still  being  debated  (Kerr,  1985). 

Another  feature  of  the  model,  which  has  received  little 
attention  is  the  collision  probability,  P  .   In  the  literature 


For  example,  Carlson  (1980)  used  S  =  0.22  for  collisions 
involving  H  atoms,  based  on  experiments  (Brackmann  and  Fite, 
19  61)  in  which  a  H  beam  impinged  upon  a  liquid-nitrogen- 
cooled  copper  plate.   The  beam  was  so  dense  that  recombination 
of  atoms  to  produce  H2  was  a  significant  source  of  loss  of 
the  atoms  from  the  beam.   This  experiment  does  not  appear 
to  closely  approximate  the  physical  system. 
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P   is  consistently  expressed  as  a  function  of  the  optical 
thickness,  t,  the  definition  of  which  is  discussed  in  Section 
III.B.   Dennefeld  set  P   equal  to  t,  while  recent  works 
(Carlson,  1980;  Ip,  1983)  use 


P   =   1  -  e"T 
c 


This  seems  to  be  a  simplification  of  the  collision  process, 
since  one  would  expect  P   to  depend  on  such  trajectory-related 
factors  as  velocity  of  the  molecule  relative  to  the  ring 
particles  at  the  location  of  passage,  and  the  path  length 
through  the  ring  (Durisen  (1984)  states  that  P   will  depend 
on  optical  thickness,  and  the  "slant  path  of  an  ejectum 
through  the  rings,"  but  does  not  elaborate).   Also,  the  above 
expressions  for  P   do  not  account  for  multiple  collisions 
during  one  passage  through  the  rings  by  a  molecule.    The 
existence  of  multiple  collisions,  if  significantly  probable, 
has  implications  both  for  density  of  the  atmosphere  (there 
are  more  "sticking"  opportunities) ,  and  for  its  spatial 
distribution. 

The  aim  of  this  paper  is  to  elevate  the  theory  of  colli- 
sions between  ring  particles  and  atmospheric  molecules  to  a 
level  of  sophistication  on  par  with  that  of  ballistic 


Note  that  nowhere  in  Dennefeld 's  model  are  multiple 
collisions  considered.   The  impact  of  this  omission  is 
discussed  later. 
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transport  models  now  utilized  in  the  study  of  ring  atmos- 
pheres.  A  generalized  model  for  determining  P   is  developed 
in  Section  II,  while  in  Section  III  its  applicability  to 
the  actual  ring  environment  is  discussed.   Section  IV  is  a 
mathematical  analysis  of  the  expression  derived  for  P  . 
Finally,  in  Section  V  the  theory  is  applied  to  the  special 
case  of  low  energy,  isotropic  production  processes  (e.g., 
sublimation)  to  illustrate  inferences  that  can  be  made  on 
the  spatial  extent  of  the  ring  atmosphere,  based  on  the  fre- 
quency of  collisions. 
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II.   THEORY 

In  this  section,  the  general  problem  of  a  molecule  passing 
through  a  translating  target  consisting  of  randomly  dis- 
tributed macroscopic  particles  is  addressed.   The  proba- 
bility that  it  collides  with  at  least  one  particle  is 
calculated . 

Consider  a  region  in  space  of  large  extent,  L,  in  two 

dimensions,  and  of  a  smaller  thickness,  h.   This  slab  is 

populated  by  noninteracting  spherical  particles  which  are 

randomly  distributed  (spatially)  and  fixed  with  respect 

to  the  slab.   Assume  that  the  size  of  the  spheres  is  limited 

to  radii  in  the  interval  [r    ,  r     <<  h]  upon  which  the 

min   max        L 

radii  may  be  distributed  either  discretely  or  continuously. 

The  number  density  of  particles  of  size  r  is  denoted  by 

n(r).    Finally,  the  slab  as  a  whole  translates  in  a  direction 

normal  to  its  thickness  with  constant  velocity,  v  . 

1         P 

Now  suppose  that  a  molecule  traveling  along  a  straight 
path  with  speed  v   passes  through  the  slab.   Assume  that  the 
molecule  is  much  smaller  than  the  smallest  particle  in  the 
slab,  and  that  it  only  interacts  with  the  particles  through 
mechanical  collision.   The  relationship  between  the  molecule 


7  -3 

For  discrete  distributions  n(r)  has  units  m   .   For 

continuous  distributions  it  has  units  m~4,  so  that  n(r)dr 

is  interpreted  as  the  number  of  particles  per  itW  in  a  radius 

increment  [r,r+dr] . 
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and  the  slab  can  be  conveniently  described  using  direction 
cosines  as  shown  in  Figure  2(a).   Notice  that  the  x-axis  is 
chosen  to  be  parallel  to  v  . 


■■>   Y 


h 


m     m   p 


(a)   v  ft   0 


(b)  v'p  =    0 


Figure  2.   Alternate  Frames  of  Reference 

Calculation  of  the  probability  of  the  molecule  colliding 
with  at  least  one  particle  as  it  passes  through  the  slab, 
henceforth  denoted  as  P  ,  can  be  accomplished  in  the  frame 
of  reference  shown  in  Figure  2(a).   This  procedure,  however, 
is  mathematically  tedious  (see  Appendix  B) .   The  problem  is 
simplified  by  converting  to  a  frame  of  reference  in  which 
v   equals  zero  (Figure  2  (b) )  .   This  approach  is  detailed 
below. 

Suppose  that  the  conversion  to  the  frame  of  reference 


with  v   equal  to  zero  has  already  been  accomplished  (see 


Appendix  C  for  details).   For  simplicity,  let  the  path  of 
the  molecule  be  normal  to  the  slab,  i.e.,  y1  equals  zero,  and 
n(r)  be  discrete  (extension  of  the  results  for  this  case  will 
be  generalized  to  cases  with  continuous  r  and  arbitrary 
molecular  path) .   The  molecule  will  collide  with  any  parti- 
cle of  radius  r.  which  has  its  center  anywhere  in  the  cylin- 
drical volume  (collision  region)  whose  axis  coincides  with 

g 

the  molecule's  path  as  shown  in  Figure  3.    The  probability 


h 


molecular  path  (normal  to 
the  slab) 


collision  region 


Figure  3.   Collision  Region  for  Particles  of  Radius  r 


that  a  collision  does  not  occur  with  particles  of  size  r. 

£  1 

is  equivalent  to  the  probability  that  no  particles  of  that 


Note  that  it  was  tacitly  assumed  that  the  molecule  is 
essentially  a  point  mass.   This  enables  one  to  define  the 
collision  region  as  was  done  above.   An  equivalent  way  to 
conceptualize  the  collision  cylinder  is  to  let  the  mole- 
cule have  radius  r,  and  the  particles  be  point  masses.   The 
molecule  then  collides  with  any  particles  it  "sweeps  out" 
along  its  path.   Of  course  the  resulting  collision  region 
again  is  a  cylinder  of  radius  r. 
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size  have  centers  in  the  collision  region.   Denote  this 
probability  as  pno(r-)«   Since  the  particles  were  assumed  to 

not  interact,  the  P   (r.)  for  all  i  are  mutually  independent. 

no   1  J  c 

This  allows  one  to  write  for  discrete  r: 


P    =   1  -  n  P   (r.),   over  all  r.  .  (II. 1) 

i 


The  problem  is  to  find  an  expression  for  P   (r.) . 
r  r  no   l 

Under  the  assumptions  of  the  model,  the  uniform  proba- 
bility distribution  governs  the  location  of  each  particle, 
i.e.,  the  probability  that  a  particle  will  be  found  in  some 
portion  of  the  slab  is  proportional  to  the  volume  of  that 
portion.   For  example,  consider  one  particle  of  size  r.. 
The  probability  that  it  is  not  in  the  collision  region  is 


2      2 

L  h  -  iTr  .  h 

Pno(l)         L2h  1 


where  the  numerator  is  volume  of  the  slab  minus  the  collision 
region,  and  the  denominator  is  the  volume  of  the  slab.   Since 
the  location  of  each  particle  in  the  slab  is  independent  of 
the  location  of  other  particles,  the  probability  that  any 
two  specific  particles  of  size  r.  are  not  in  the  collision 
region  is 


2      2 

L  h  -  TTr  ■  h  „ 

Pno(2)      (~"7V 
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Extending  this  line  of  reasoning,  it  is  clear  that 


2      2 

L  h  -Tir  h  N(r  .  ) 

Pnn(r.)   =   ( 5 — )    X   , 

n°   X  L2h 


where : 


N(r.)   =   n(r.)L2h 


is  the  total  number  of  particles  with  radius  r.  in  the  slab. 


P   (r.)  can  be  rewritten  as 
no   i 


2 

TTr  .  n  (r  .  )  h  N  (r  .  ) 

Pno(ri>   =   (1  -   N(r.)   '    "  (II-2) 


If  N(r.)  is  large,  then 


P   (r.)   =   exp{-iTr2n(r .  )h}  .  (II. 3 

no   i        ^ii 


Substituting  into  equation  (II. 1)  with  N(r.)  large  for  all 

r.  yields 

l  2 


2 

P       1  -  n  exp{-7Tr  .n  (r  .  )h} 

c  ^11 

l 


or 


P    =   1  -  exp{-7Th  Jr2n(r.)}  ,  (II. 4) 

C  c  h      l     l 

1 
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recalling  that  r  is  discretely  distributed,  and  the  molecular 
path  is  normal. 

Extending  this  result  to  the  case  when  the  molecule's 
path  makes  an  angle  y '  with  a  normal  from  the  slab  is  easily 
accomplished.   Notice  that  the  only  difference  between  this 
case  and  that  of  the  normal  path  is  that  the  length  of  the 
molecular  path  (hence  the  length  of  the  collision  cylinder) 
increases  from  h  to  h/ | cos  y'I-   Therefore, 

Pc   ■   1  -  exPflcosVi  lri"'ri»>  <II-5> 

I       I  1 

with -r  discrete  and  the  molecular  path,  not  necessarily  normal 

Extension  to  the  case  when  r  is  continuous  is  affected 
by  replacing  the  summation  in  equation  (II. 5)  with  integra- 
tion (see  the  first  footnote  on  page  9  to  see  the  new  inter- 
pretation required  for  n(r)).   This  yields 


r 
,       max  ? 

P    =   1  -  exp{n ,  i  i  r  n(r)dr}  ,  (II. 6) 

c  r   cos  Y 

i     r  i   r 

min 


where  r  is  continuous  and  the  molecular  path  is  not  neces- 
sarily normal  to  the  slab. 

Equations  (II.  5)  and  (II. 6)  represent  solutions  to  the 
problem  when  working  in  a  frame  of  reference  for  which  v 
equals  zero.   Often  problems  are  not  initially  formulated 
in  this  particular  frame.   Appendix  C  shows  how  to  convert 
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these  results  back  to  the  original  frame  of  reference  (shown 
in  Figure  1).   Substituting  from  equation  (C.8)  yields 


Pc  -  l   -exP{|coshylV1'2(^)COSa+'^)2  £rin(ri)}        (II-7» 


for   r   discrete,    and 


P       =1  -exp{- 


c  L     I  cos   y 


Vv  v      ?        max      o 

l-2(^£)cosa+(^£)^    /  r    n(r)dr} 

m  m  r    . 


m  m  r 

mm 


(II. 8 


for  r  continuous.   Equations  (II. 7)  and  (II.  8)  are  more  con- 
venient than  equations  (II. 5)  and  (II. 6)  for  some  applications 
It  is  interesting  to  note  that  while  the  length  of  the 
molecular  path  in  the  frame  of  reference  with  v   not  equal 
to  zero  is  h/|cos  y|,  the  "effective  path  length,"  pertinent 
to  the  collision  problem,  is 


,        /       V  V    9 

L  --   =   i r\/l  -2(-£)cos  a   +  (-£)   ,  (II. 9) 

ef f      cos  Y  V      v  v 

1     '  '  V       m  m 


which  accounts  for  the  relative  motion  of  the  molecule  versus 
the  ring  particles. 

Equations  ( II . 5) -  ( II . 8)  represent  the  fruition  of  this 

derivation  of  P  .   Table  1  gives  a  review  of  the  tacit  and 

c  a 


explicit  assumptions  of  the  model 
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TABLE  1 
ASSUMPTIONS  OF  THE  THEORETICAL  MODEL 


Assumptions 


Comments 


spherical  particles 

random  spatial  distribu- 
tions of  the  particles 

h  <  L 


r  .   <.  r  <  r 

mm  —   —  max 


r     <<  h 

max 


r  .   >>  r   ,    , 
mm  ;-.    molecule 


molecule  enters  and  leaves 
slab  through  its  large 
(area  of  Ir-)    faces 

slab  translates  with 
velocity,  v  ,  normal  to 
its  thickness,  h 

-»-      ->• 

v   and  v   constant 
p       m 

particles  fixed  relative 
to  the  slab 

number  of  particles  in  the 
slab  is  large 


mathematically  convenient 

justifies  using  the  uniform 
probability  density  function 

facilitates  application  of 
the  model  to  the  ring-s 

no  significance  beyond 
"quantifying"  particle  size 

effectively  excludes  colli- 
sions outside  the  nominal 
confines  of  the  slab 

allows  the  molecule  to  be 
treated  as  a  point  particle 

facilitates  application  of 
the  theory  to  the  rings 


facilitates  application  of 
the  theory  to  the  rings, 


mathematically  convenient 
mathematically  convenient 


justifies  exponential 
approximation  in  probability 
calculation 


Is  the  theoretical  model  utilized  here  relevant  to  condi- 
tions in  Saturn's  rings?   This  question  is  discussed  in  the 
next  section. 
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III.   APPLICATION  TO  SATURN'S  RINGS 

A.   APPLICABILITY  OF  THE  THEORETICAL  MODEL 
1 .   General  Comments 

Before  proceeding  with  analysis  of  collisions  between 
ring  particles  and  atmospheric  molecules,  the  validity  of  the 
generalized  "slab"  model  of  Section  II  must  be  established 

vis-a-vis  the  actual  ring  environment.   The  essential  assump- 

9 
tions   of  the  model  are: 

(1)  the  particles  in  the  slab  are  spherical. 

(2>  the  radius  (size)  of  any  particle  in  the  slab  is 
in  an  interval  [rmj_n,  rmax]  where  rmj_n  is  much 
greater  than  the  size  of  the  molecule  and  rmax 
is  much  smaller  than  the  slab  thickness,  h. 

(3)  the  spatial  distribution  of  particles  within  the 
slab  is  random,  i.e.,  the  location  of  each  particle 

is  governed  by  the  uniform  probability  density  function 
with  all  particles  mutually  independent  (in  a  sta- 
tistical sense)  . 

(4)  the  number  of  particles  within  the  slab  is  large 
enough  for  the  exponential  approximation  (equation 

(II. 3) )  to  be  valid. 

(5)  the  velocity  of  the  slab,  v  ,  is  constant. 

i  p 

(6)  the  velocity  of  the  molecule,  v  ,  is  constant. 

(7)  the  slab  particles  are  fixed  relative  to  the  motion 
of  the  slab  as  a  whole. 


9 

Some  of  the  assumptions  listed  in  Table  1  were  adopted 

to  facilitate  application  of  the  model  to  the  rings.   These 
assumptions  have  no  physical  significance,  so  they  are  not 
discussed  here. 
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If  a  "physical  slab"  can  be  defined  in  the  rings  which 
reasonably  satisfies  the  above  conditions,  then  the  appli- 
cability of  the  model  is  verified.   Note,  however,  that  the 
physical  slab  is  an  abstract  concept,  since  it  is  related  to 
the  passage  of  a  specific  molecule,  on  a  specific  trajectory. 

To  facilitate  the  discussion  of  the  above  seven 
assumptions,  they  can  be  divided  into  three  categories,  those 
related  to: 

(1)  particle  properties  (assumptions  (1)  and  (2)). 

(2)  particle  distribution  properties  (assumptions  (3) 
and  (4)). 

(3)  velocity-related  properties  (assumptions  (5)-(7)). 
Below  each  of  these  categories  is  discussed  with  respect  to 
Saturn's  A  and  B  rings.   The  other  rings  will  not  be  con- 
sidered, since  they  are  not  expected  to  contribute  much  to 

the  ring  atmosphere  (owing  to  their  low  masses) .   Observational 
data  referenced  in  Section  I  supports  this  exclusion. 
2 .   Particle  Properties 

Individual  ring  particles  have  not  yet  been  observed, 
thus  absolute  verification  of  the  assumption  of  spherical 
shape  cannot  be  accomplished.   However,  such  a  shape  is 
adopted  ubiquitously  in  the  development  of  ring  models. 
Observational  data  does  not  suggest  that  this  traditional 
approach  is  invalid  (Weidenschilling  et  al.,  1984),  in  fact 
the  Voyager  radio  occultation  experiment  (Tyler  et  al . ,  1983) 
eliminates  the  possibility  of  very  elongated  particles. 
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Therefore,  at  present  the  most  reasonable  shape  to  attribute 
to  ring  particles  is  spherical. 

Most  ring  particles  have  sizes  in  a  continuous  radius 
interval  (Cuzzi  et  al .  ,  1984;  Esposito  et  al.,  1984; 
Weidenschilling  et  al.,  1984) 

1  cm   <   r   <   5m. 

Clearly,  r     is  much  greater  than  the  size  of  a  H~0  or  H 
mm  2. 

molecule.   The  thickness,  h,  of  the  A  and  B  rings  is  nominally 

not  larger  than  200  m,  but  locally  may  be  as  small  as  10  m 

(Cuzzi  et  al.,  1984).   The  requirement  that  r     be  much  less 

^  max 

than  h  should  not  create  a  problem  for  application  of  the 
collision  theory. 

3 .   Spatial  Distribution  of  the  Ring  Particles 

The  simplest  way  to  describe  the  rings  is  as  a 
collection  of  mutually  colliding  particles  whose  trajectories 
are  dominated  by  the  gravitational  influence  of  Saturn.   Under 
such  conditions  the  notion  of  a  completely  random  particle 
distribution  throughout  the  rings  is  reasonable.   The  physi- 
cal slab  could  then  be  defined  as  any  translating  region 
containing  the  molecular  path,  provided  that  the  velocity- 
related  assumptions  are  valid  for  that  specific  region. 
Choose  the  largest  such  region  consistent  with  the  other 
assumptions  to  be  the  physical  slab.   As  will  be  shown  in 
the  discussion  of  velocity-related  requirements,  this  region 
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will  have  a  considerable  volume.   This  is  desirable,  for  even 
if  n(r)  is  quite  small  (as  limited  measurements  have  indi- 
cated (Marouf  et  al . ,  1983)),  the  requirement  that  the  number 
of  particles  in  the  slab  be  large  is  likely  to  be  satisfied. 

In  actuality  the  rings  exhibit  considerable  fine 
structure  (i.e.,  variations  in  optical  thickness)  on  radial 
scales  down  to  -1  km  (Esposito  et  al . ,  1984).   Apparently 
the  primary  cause  of  this  fine  structure  is  linked  to  the 
gravitational  effects  of  other  bodies  in  the  Saturnian  system, 
especially  the  moons.   Gravitational  resonances  are  respon- 
sible for  the  creation  of  zones  where  stable  particle  orbits 
are  precluded  (e.g.,  the  Cassini  Division),  and  may  generate 
density  waves  within  the  rings.   The  net  effect  of  fine  struc- 
ture on  the  collision  theory  is  to  limit  the  scale  over  which 
random  spatial  particle  distribution  is  correct.   This  can  be 
best  understood  by  considering  the  following  examples. 

In  the  first  example,  the  slab  is  divided  into  two 
regions  as  shown  in  Figure  4.   Suppose  that  the  particles  in 


molecular        Y 
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Note :   For  simplicity, 
consider  the  frame  of 
reference  with  v   =  0, 
and  require  the  molecular 
path  to  be  parallel  to 
a  slab  edge. 


Figure  4.   Model  for  Analyzing  the  Effect  of  Nonrandom 
Particle  Distribution 
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each  region  are  constrained  to  remain  within  that  region; 
this  represents  the  departure  from  truly  random  spatial  dis- 
tribution of  the  particles.   Let  n, (r)  equal  n~(r)  with 
r  discrete.   What:  is  the  probability  that  the  molecule 
suffers  at  least  one  collision  when  passing  through  the  entire 
slab?   Attack  this  problem  by  reformulating  it  as  passage 
through  two  smaller  slabs  shown  in  Figure  5  as  shaded  regions. 


b 


Figure  5.   The  Two  Slabs  used  in  Reformulating  the  Problem 


From  Section  II,  the  probability  of  no  collision  with  parti- 
cles of  size  r.  for  the  slab  of  thickness  b  is 

1 


Pno(b)(ri)   = 


1  - 


2  /   x  ,     b 

nr.n  r.   -, r 

l    l    cos  y 


N(ri) 


N(r 


where: 


N(r.)   =   n(r.)L  b 
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For   all    r .    this   becomes 


no(b) 


n    (1 
i 


2    /       x  ,  b 

irr  .  n  (r  .  )  (-1 

1        1        cos   y 


N(ri) 


N(r  .  ) 

1 


which    for    large   N(r.)    becomes 


no(b) 


=      II    exp{-7rr  .  n  (r  .  ) 
L  11 

1 


cos   y 


7-r)} 


r  — TTb  r  2       ,  .   i 

exp  n j-i-  )    r.n(r.)} 

c      cos   y         •       11 


III.l) 


Similarly,  for  the  slab  of  thickness  (h-b) 


no(h-b) 


exp{- 


TT(h-b)   r   2 


cos  y 


I    rfn(r.) }  (III. 2) 

r       1    1 


for  n(r.)L  (h-b)  large.   The  probability  of  collision  for  the 
entire  passage  through  the  original  slab  is 


Pc      -1    Pno(b)  Pno(h-b)  * 


(III. 3) 


Substitution  yields 


„  !  r       _TTb       r     2    ,      .  ■.         r   — tt  (h-b)    r     2    .     .  1 

P       =     1  -  exp{-i rr  )   r.n(r . )  }exp\-i s A-  )   r.n(r.)}      , 

c  r     cos  y        ■      11  cos  Y        •      11 
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or 


P    =   1  -  exp{  I   ~  ",  I  I  r^n(r.)}  .  (III. 4 

c  c      cos  Y    .   1    i 


l 

This  result  is  the  same  as  that  obtained  when  a  random  spatial 
distribution  over  the  entire  slab  is  assumed.   Careful  con- 
sideration reveals  an  important  qualification  to  this  conclu- 
sion.  Notice  that  the  assumption  of  exponential  form  is  less 

tenable  in  the  "nonrandom"  case,  since  there  it  is  required 

2  2 

that  both  n(r.)L  b  and  n(r.)L  (h-b)  be  large,  while  in  the 

2 
completely  random  case  n(r.)L  h  must  be  large.   Extending  this 

line  .of  reasoning  to  cases  where  the  original  slab  is  further 

subdivided  into  regions  to  which  the  particles  are  confined, 

it  is  clear  that  at  some  point  the  exponential  approximation 

will  certainly  fail. 

In  the  second  example,  the  slab  is  again  divided  into 

two  regions  (Figure  4) ,  but  now  n, (r)  is  not  equal  to  n_(r) . 

Following  the  approach  taken  above,  the  collision  probability 

is 


P   =  l-exp{,  "ub,,  I  r2n,(r.)}exp{r1T{h-^\    Y  r2n0(r.)}  , 
c        r  cos  Y   ■   ill       cos  Y    •   1  2  T 


or 


P^  =  1  -exp{  i  ~J    ,  |  [b  T  r2ni(r.) 

C  c        COS  Y        •    1  1    1 


cos  y 

1 


+  (h-b)  I    r2n0(r .)] }  ,  (III. 5 

j  J    2  2 
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2  2 

where  both  n, (r.)L  b  and  n2(r.)L  (h-b)  must  be  large  enough 

to  justify  the  exponential  approximation.   The  collision 

theory  is  here  valid  in  principle,  but  its  application  is 

more  difficult. 

From  this  point  on,  the  A  and  B  rings  will  be  assumed 
to  satisfy  the  spatial  distribution  conditions.   However,  it 
must  be  kept  in  mind  that  fine  structure  effects  may  inter- 
ject inaccuracies  into  calculations. 
4 .   Velocity-related  Assumptions 

Reviewing,  the  three  velocity-related  assumptions 

are:   v   is  constant;  v   is  constant;  and  the  particles  in  the 
p  m  '  t- 

slab  -are  fixed  relative  to  v  .   Each  in  turn  is  addressed 

P 

below,  but  the  validity  of  the  concept  of  slab  velocity  when 
applied  to  the  ring  system  must  first  be  demonstrated.   As 
discussed  in  Section  III. A. 3,  choosing  the  physical  slab  to 
be  as  large  as  possible  is  advantageous.   By  means  of  a  simple 
example,  it  can  be  shown  that  a  single  velocity  well  approxi- 
mates the  orbital  motion  of  all  of  the  particles  in  even  a 
quite  large  slab  (this  velocity  corresponds  to  v  ) .   Suppose 
that  a  slab  with  L  equal  to  10,000  km  is  centered  in  the  B 
ring  as  shown  in  Figure  6.   If  v   is  set  equal  to  the  parti- 
cle velocity  at  C,  v  ,  how  much  error  is  inherent  to  the 
system?   First,  consider  the  magnitude  of  v  .   The  location  in 
the  slab  with  the  greatest  deviation  in  speed  relative  to  v 
is  either  point  A  or  B.   In  general  the  orbital  speed  for  a 
ring  particle  in  the  Saturnian  system  is 
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R   =  99,790  km 

A 


Saturn ' s 


Center 


R  =92,000  km 
m 


R_  =  109,79  0  km 
B 


R 


RAT  =  117,580  km 
N 


B  RING 


Figure  6.   A  Slab  with  L  =  10,000  km  Centered  in  the 
B  Ring 


v   =   /GM  /R   =   iili   [kms  1]     .  (III. 6 


/R 


Applying  this  equation  at  A,  B,  and  C  yields 


and 


vA   =   19.496  kms  X  , 


vB   =   18.5  87  kms  1  , 


v    =   19.025  kms 
c 


-1 


Thus,  all  of  the  particles  in  the  slab  orbit  Saturn  with  a 
speed  in  the  range  [18.587  kms   ,  19.496  kms   ] ;  the  greatest 
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deviation  in  assigning  v   (equals  v  )  to  all  of  the  particles 

c  p 

is  2.5%.   Now,  consider  the  direction  of  v  .   The  point  of 
greatest  deviation  in  direction,  J> ,  occurs  at  point  E  (the 
angle  of  deviation  at  point  F  is  equal  to  that  at  E) .   Figure 
7  gives  the  computation  of  $;  the  results  given  there  demon- 
strate  that  the  angular  error  in  assigning  v   to  all  of  the 


of 
Saturn 


.   -1,  5000  , 
*  =  tan   (99,790} 

-2 
<b   =  5.006  x  io    rad 


R_  =  9  9,915  km 

E 

v„ I  =19.4  85  kms 
E 


-1 


vn  =0.97  5  kms 
Ex 


-1 


v^  =19.4  60  kms 
Ey 


-1 


Figure  7.   Calculation  of  the  Maximum  Angle  of  Deviation,  $ 


particles  is  quite  small.   The  above  example  verifies  the 
concept  of  a  slab  velocity,  v  ,  for  large  slabs  in  the  B 
ring  (similar  results  also  hold  for  the  A  ring) .   The  concept 
applies  even  better  for  smaller  slabs. 

Accepting  that  v   is  a  meaningful  concept  with  regard 
to  physical  slabs  in  the  A  and  B  rings,  the  requirement  that 
it  be  constant  can  now  be  examined.   Once  defined,  the  slab 
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will  translate  as  shown  in  Figure  8.   Upon  reflection,  the 
slab  velocity  cannot  be  strictly  constant.   Although  |v 

vH 


po 


Cente 


slab  at 

time  t,  >  t 
1    o 


slab  at 
time  t 


Figure  8.   Translation  of  the  Slab  over  Time 


equals  |v  , | ,  the  direction  of  the  velocity  vector  changes 
with  time.   The  angle  cj)  which  quantifies  the  change  in  direc- 
tion is  given  by 


v 


$   = 


(III. 7 


For  a  slab  centered  in  the  B  ring  (as  in  Figure  6) , 


<J>   =   1.8156  xl0"4t  , 


(III. 8) 
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where  t  is  in  seconds,  <j>  in  radians.   The  time  required  for 

-2 

a  variation  of  1.7453  xlO    rad  (1°)  is  96.1  s.   This  example 

suggests  that  for  the  A  and  B  rings,  the  assumption  of  con- 
stant  v   is  approximately  correct,  if  the  time  of  passage 

ir 

for  a  molecule  through  the  ring  is  short.   The  requirement 
will  be  met  for  most  passages  (i.e.,  over  most  y) ,    since  the 
thickness  of  the  ring  is  so  small.   To  see  this,  then  sup- 
pose that: 

(1)  h  =  100  m, 

(2)  R  =  105,000  km  (molecule  passes  near  the  center  of 
the  B  ring) , 

(3)  I v  I  =20  kms"1.10 
..   '  m ' 

Over  what  range  of  y    (recall  Figure  2)  will  cj>  be  less  than  1°; 
or  phrased  differently,  over  what  range  of  y  will  the  time 
of  passage  of  the  molecule  be  less  than  96  s?   In  general  the 
time  for  passage  for  the  molecule  is 


t    =    r-£ r   .  (III.  9) 

vjcos  yI 


Substituting  the  given  information,  requiring  t  to  be  equal 
to  96  s,  and  solving  for  y  yields 


96   =   10° 


20x10  | cos  y 


The  lowest  kinetic  energy  of  a  stable  orbit  (measured 
in  a  frame  of  reference  fixed  with  respect  to  Saturn)  for  an 
atmospheric  particle  which  originates  in  the  A  and  B  rings, 
as  it  passes  through  those  rings,  is  -20  kms--'-. 
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cos  y I   =   5.21  x 10  3  , 


y   =   89.997°,  90.003°  . 
Therefore,  in  this  representative  case  the  range  of  y    is 
0   £   y   <   89.997°  ,    90.003°   <  y  180°  ; 

virtually  any  y  meets  the  requirement. 

-»■ 
The  requirement  that  v   be  constant,  as  in  the  case 
^  m 

■> 
of  v  ,  cannot  be  strictly  met,  since  the  trajectory  of  any 

molecular  orbit  will  have  some  curvature.     However,  the 

extreme  thinness  of  the  rings  allows  the  condition  to  be 

approximately  satisfied,  so  this  aspect  of  the  theoretical 

model  can  be  applied  to  Saturn's  rings  with  no  problem. 

That  the  particles  on  the  slab  are  fixed  with  respect 

to  v  cannot  be  strictly  satisfied,  either.   Direct  measure- 

ments  of  the  relative  velocities  of  the  ring  particles  have 

not  been  accomplished,  but  various  indirect  methods  yield 

consistent  results:   ~1  cms    (Weidenschilling  et  al . ,  1984). 

Compared  to  the  range  of  orbital  speeds  for  particles  in  the 

A  and  B  rings 


16.6  kms     <   v    <   2  0.3  kms    , 

P   - 
(outer  A  ring)  (inner  B  ring 


Note  that  the  gravitational  interaction  between  the 
molecule  and  ring  particles  has  been  neglected,  a  reasonable_ , 
assumption  since  the  total  mass  of  the  rings  is  less  than  10 
that  of  Saturn  (Null  et  al . ,  1981). 
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calculated  using  equation  (III. 6),  the  relative  motion  of 
the  particles  is  insignificant  over  reasonable  passage  times 
(less  than  or  equal  to  -100  s,  as  discussed  above).   Thus, 
in  this  regard  the  theoretical  model  suffers  no  important 
restriction  in  applicability  to  the  rings. 
5 .   Summary 

The  above  discussion  has  demonstrated  that  the  condi- 
tions upon  which  the  theoretical  model  is  based,  are  gener- 
ally well-satisfied  in  the  ring  environment.   The  only 
potential  source  of  significant  concern  relates  to  assumption 
of  random  particle  distribution  and  the  validity  of  the  exponen- 
tial ..approximation  (Section  III. A. 3).   Unlike  the  other  assump- 
tions of  the  model,  these  cannot  be  unamibguously  shown  to 
hold  for  Saturn's  rings.   On  the  other  hand,  they  cannot  be 
shown  to  fail,  and  have  considerable  intuitive  appeal.   With 
this  caveat  noted,  the  theoretical  model  is  adopted  for  analy- 
sis of  collisions  between  ring  particles  and  atmospheric 
molecules . 

B.   THE  UTILITY  OF  OPTICAL  THICKNESS  IN  APPLICATIONS 

Having  shown  that  the  theoretical  model  can  be  viably 
applied  to  the  atmosphere  in  the  vicinity  of  the  A  and  B 
rings,  calculation  of  P   using  the  results  of  Section  II  is 
justified.   Since  particle  radius  is  continuously  distributed, 
and  effectively  confined  to  an  interval  {0.01  m,5  m]  as  dis- 
cussed earlier,  the  collision  probability  is  (from  Section 
ID 
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h      5       5 

Pc   =   1  -  exp{|cos^,|  /      r  n(r)dr}  ,  (II. 6) 


when  working  in  a  frame  of  reference  with  v   equal  to  zero, 
and  is 


,     /    v         v   p   5      „ 

P   =  l-exp{,  ~1Tn    |  Wl-2(-E)cos  a +  {-&■)  /   '   rn(r)dr}  , 

C  r    COS  Y   V       V  V        o  m 

'     '  '  v      m         m     0.01 


II. 


when  working  in  a  frame  of  reference  where  v   is  not  equal  to 
zero.   The  value  for  h  is  well-constrained,  while  the  param- 
eters v  ,  v  ,  a,  y.    and  y1  depend  on  the  location  of  a  mole- 
m   p      '  r 

cule's  passage  and  its  trajectory  (both  of  which  can  be 
determined  using  Keplerian  mechanics).   If  n(r)  is  known 
throughout  the  rings,  then  the  calculation  of  P   using  the 
above  equations  will  be  straightforward.   Although  the  func- 
tional form  of  n(r)  is  believed  to  be 


n  (r)   =   n  r  p  ,   p   >   0  , 


where  p  may  differ  over  r,  its  precise  expression  has  only 
been  determined  for  a  few  locations  in  the  A  and  C  rings 
(Marouf  et  al . ,  1983).   Fortunately,  this  obstacle  to  the 
application  of  equations  (II.  6)  and  (II.  8)  can  be  overcome 
by  using  the  concept  of  "optical  thickness." 


39 


12 
Optical  thickness   ,  t,  is  defined  in  an  operational 

manner:   if  an  electromagnetic  beam  traverses  a  layer  of 

particles  at  an  angle  y  to  the  layer  normal,  and  the  layer 

has  an  optical  thickness  of  t,  then  the  beam's  intensity  is 

—  -[■  /co  s  u 
reduced  by  a  factor  of  e         .   In  general  t  will  depend 

on:   particle  composition,  size,  and  shape;  particle  density; 
layer  thickness;  and  wavelength  of  the  incident  beam.   De- 
tailed analysis  of  this  phenomenon  is  beyond  the  scope  of  this 
paper  (see  van  de  Hulst  (1981)  and  Chandrasekar  (1960)  for 
thorough  analyses  of  light  scattering  by  media  consisting  of 
small  particles).   However, the  essential  attributes  of  a  com- 
plete- analysis  of  t,  can  be  concisely  summaried,  generally 
following  Cuzzi  et  al .  (1984). 

For  spherical  particles  in  a  layer  of  thickness,  h,  the 
optical  thickness  is  given  as 


t(A)   =   /  Q  (r,A)Trr2hn(r)dr  ,  (III. 10) 


where : 


wavelength  of  incident  radiation, 


and 


12 

Optical  thickness  is  often  referred  to  as  "optical 

depth."   The  terms  are  used  interchangeably  in  literature 
concerning  planetary  rings.  Also,  since  t  depends  on  wave- 
length, one  must  be  careful  when  comparing  the  results  of 
different  studies. 
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Q  (r,A)   =   extinction  cross  section  (a  measure 
of  the  portion  of  incident  energy 
removed  from  the  beam  by  absorption 
and  scattering) . 

When  r  is  much  greater  than  A ,  then 

!1 ,  for  incoherent  incident  radiation. 
2,  for  coherent  incident  radiation. 

Considering  visible  wavelengths  with  respect  to  the  predominant 
particle  sizes  in  the  A  and  B  rings,  i.e.,  [0.01  m,5  m] ,  the 
above  values  for  Q   apply.   Therefore,  for  the  A  and  B 
"rings  the  visible  optical  thickness  when  measured  with  an 
incoherent  signal  (e.g.,  stellar  occultation)  is 


x   =   /  Trr2hn(r)dr  .  (III. 11) 


Equations  (II. 6)  and  (II.  8)  can  be  rewritten  as 


P    =   1  -  e-T/'COS  Y'l  (III. 12 

c 


Vv  V   ? 

l-2(-£)cos  a  +  (-£)  }  (III. 13 

V  V 
m  m 


Of  course,  if  optical  thickness  measurements  based  on  coherent, 
visible  signals  are  used  to  compute  x,  the  resulting  expres- 
sions for  P   will  differ  from  equations  (III. 12)  and  (III. 13) 
by  a  factor  of  0.5  in  the  exponent  of  e. 
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Measurements  of  i  for  Saturn's  rings  abound.   To  apply 
equations  (III. 12)  and  (III. 13)  to  the  rings,  one  must  only 
be  cognizant  of  the  importance  of  the  wavelength  at  which 
t  was  determined.   Keeping  this  in  mind,  analysis  of  P   can 
proceed. 
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IV.   GENERAL  ANALYSIS  OF  THE  EXPRESSION  FOR 
THE  COLLISION  PROBABILITY,  Pn 

A.   RESTRICTIONS  ON  THE  PARAMETERS 

When  studying  Saturn's  rings,  the  frame  of  reference 
usually  adopted  features  Saturn  as  a  fixed  central  body 
orbited  by  ring  particles  following  (circular)  Keplerian 
trajectories.   It  is  also  the  reference  frame  of  most  general 
interest  when  examining  the  ring  atmosphere,  since  ballistic 
transport  models  are  based  on  basic  orbital  mechanics.   For 

this  reason  only  equation  (III. 13)  is  analyzed  in  this 

13 

sect-ion. 

Recall  that  the  collision  probability  is  given  as 


Vv  v   ? 

l-2(-£)cos  a  +  (-£)  }  (III. 13) 

vm         vm 


Clearly  P   is  a  function  of  parameters  related  to  the  mole- 
cule's orbit  (v  ,a,y),  and  those  that  characterize  the 

m 

translating  slab  (x,v  ).   Note  that  a    and  y  are  not  indepen- 


dent, but  must  satisfy 


2        2 
cos  a  +  cos  y   <   1  t 


13 

That  is  not  to  say  that  the  reference  frame  with  slab 

particles  at  rest  is  not  useful  in  some  practical  applications 
(as  is  shown  in  Section  V) .   Rather,  although  the  collision 
problem  may  be  solved  in  either  frame  of  reference,  the  prob- 
lem of  calculating  molecular  trajectories  is  ultimately  ap- 
proached in  the  frame  with  v   not  equal  to  zero,  requiring 
conversion  to  that  frame.    p 
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since  they  are  direction  cosines.   An  alternative  coordinate 
system  which  removes  this  interdependence  specifies  a  mole- 
cule's path  in  terms  of  y  and  0,  henceforth  referred  to  as 
the  "orbital  inclination"  and  "orbital  azimuth."   9  is  here 
defined  as  the  angle  between  the  projection  of  v   in  the 
x-y  plane  (the  ring  plane),  and  the  x-axis  (i.e.,  the  direction 
in  which  the  slab  translates) .   Figure  9  shows  the  relation- 
ship between  0  and  the  direction  cosines,  the  mathematical 
expression  of  which  is  given  in  Appendix  D.   Employing  the 


-/  y 


Note 


-7T   <   0   <   TT 


Figure  9.   Definition  of  the  "Orbital  Azimuth,"  0 

results  of  Appendix  D,  the  collision  probability  can  be 
rewritten  as 


-T 


P   =  l-exp{- 
c         r  | cos  y 


Vv  v   ? 

l-2(-£)sin  y  cos  9  +  (-£)  } 
vm  Vm 


(IV. 1) 


Equation  (IV. 1)  proves  more  useful  than  equation  (III. 13 
in  the  subsequent  analysis. 
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Before  examining  the  expression  for  P   in  detail,  physi- 
cal considerations  which  limit  the  range  of  values  for  the 
parameters  are  outlined  below.   Following  this,  attention  can 
then  be  turned  to  parameter  values  of  most  practical  interest 
in  the  study  of  the  atmosphere  of  Saturn's  rings. 

The  optical  thickness,  t,  of  Saturn's  rings  has  been 
studied  extensively.   Recall  from  Section  III.B  that  the 
measured  value  of  t  is  dependent  on  several  variables,  most 
importantly  wavelength  and  coherence  of  the  incident  radiation, 
and  the  nominal  size  of  the  ring  particles  at  the  location 
where  the  electromagnetic  beam  is  incident  upon  the  ring  plane. 
As  shown  in  Section  III.B,  measurements  based  on  incoherent, 
short  wavelength  (e.g.,  visible)  radiation  yield  values  of 
t  that  can  be  most  easily  used  in  collision  theory  (see 
equations  (III. 11)  and  (111.13)).   A  stellar  occultation 
experiment  onboard  Voyager  2  provided  results  that  meet  the 
above  requirements  (Esposito  et  al.,  1984).   Optical  thick- 
ness was  found  to  fluctuate  radially,  varying  from: 

(1)  -0.1  to  -0.4  in  the  C  ring. 

(2)  -0.7  to  -2.5  in  the  B  ring. 

(3)  -0.1  to  -0.2  in  the  Cassini  Division. 

(4)  -0.4  to  -1.0  in  the  A  ring. 

Average  x  over  the  A-C  rings  is  given  in  Appendix  A.   Sum- 
marizing this  information,  the  nominal  optical  thickness  is: 

(1)  0.1  for  the  C  ring  and  the  Cassini  Division. 

(2)  1.5  for  the  B  ring. 

(3)  0.5  for  the  A  ring. 
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Subsequently,  in  further  analysis  of  Saturn's  ring  atmosphere, 
t  will  be  limited  to  the  interval 

0.1  <_      t      2.0  . 

The  slab  speed  (which  is  essentially  the  ring  particle 

speed  in  the  vicinity  of  the  molecule's  path  through  the 

ring)  and  the  molecular  speed  appear  as  a  ratio,  v  /v  ,  in 

the  expression  for  P  .   Call  v  /v   the  "speed  ratio."   Esti- 
c  c         p  m       r 

mates  of  H„0  molecule  ejection  speeds,  relative  to  the  source 
particle,  for  several  production  processes  are: 

(1)   ~0.3  kms    for  sublimation  (Dennefeld,  1974). 

(2>   ~1  kms    for  meteroid  bombardment  (Dennefeld,  1974). 


Estimates  for  H  atom  ejection  speed  are: 

<  3  kir 
1983) . 


(1)   <  3  kms    for  meteroid  bombardment  (Morfill  et  al., 


(2)   <  10  kms    for  photodissociation  of  H-0  (Durisen,  1984) 
Based  on  the  above  values,  estimates  for  the  speed  ratio 
range  for  molecules  newly  ejected  from  particles  in  the  A 
and  B  rings  are 


v 

0.95  <      S.      <      1.06 
—  v    — 
m 


for  H~0  molecules,  and 


v 

0.40   <   -2-   <   1.49 
—  v    — 
m 
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for  H  atoms,  noting  that  v   is  the  sum  of  v   and  the  ejection 

^       m  P 

velocity,  v  ..   The  eccentricity  of  the  subsequent  orbit, 
and  the  effect  of  collisions  complicate  precise  specification 
of  general  limits  on  the  speed  ratio,  but  the  above  should 
suffice  as  guidelines.   Therefore,  further  analysis  of  P   in 
this  paper  will  be  limited  to  a  range 


v 

0   <   -2-   <   2  . 
—  v    — 
m 


By  definition,  y    and  9  are  restricted  to  [0,tt]  and  [-tt,tt]  , 
respectively.   Considering  equation  (IV. 1),  it  can  be  seen 
that  - 


P  (9)   =   P  (-9)  , 
c         c 


and 


PC(Y)   =   Pc(tt-y)  +-*    P  (jir-y)   =   Pc(jTT+y)  . 


Noting  these  symmetries,  no  information  is  lost  by  limiting 

subseauent  analysis  of  P   to 

1  c 


0    <    9    <    TT   , 


and 


0  <_      y       <_      -   tt 
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Though  P   will  be  analyzed  over  the  above  intervals,  some 
physical  restrictions  on  values  for  orbital  inclination  and 
azimuth  are  discussed  in  the  following  paragraphs. 

When  a  given  production  process  is  considered,  the  initial 
limits  on  y    can  easily  be  determined  as  shown  in  Figure  10. 
For  example,  the  y  ""range  for  H-O  is: 


v   =  v   +  v 
m    p     ej 


z 


•+    X 


locus  of  possible 
ejection  direction 


v 


m 


1 

nn    ' 


x  v. 


^< 


e: 


■*  x 


v 


min 


v 


-1   ei 
Note:  YmA„    =  cos   (~^-J-) 

!b) 


Figure  10.   (a)   General  Relationship  Between  vej ,  vm, 

and  Vp . 
(b)   Case  for  Minimum  Value  of  y 


(1)  89°   <   y   ~   91°    for  sublimation. 

(2)  86.6°   <   y   ~   93.4°    for  meteroid  impact. 
Similarly,  for  H  atoms,  y    is  initially  confined  to: 

(1)   79.6°   <   y   ~   100.4°    for  meteroid  impact. 


(2)   53.0°   <   y   ~ 


127.0°    for  photodissociation  of  H20 
on  a  particle  surface. 
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As  collisions  occur  between  molecules  and  ring  particles, 
the  range  on  y  will  change.   If  a  collision  is  assumed  to  be 
inelastic,  then  upon  impact  the  molecule  will  rebound  with 
reduced  kinetic  energy  relative  to  the  particle.   The  overall 
effect  is  that  y  should  go  to  -jtt  as  more  collisions  occur. 

For  atmospheric  molecules  originating  on  a  ring  particle, 
the  initial  orbital  azimuth  will  certainly  be  "forward- 
directed,  "  i.e. , 

as  long  as  the  ejection  speed  is  less  than  the  particle  speed 

This  again  follows  from  the  fact  that  v   is  the  vector  sum  of 

3  m 

->■       -*■ 

v   and  v  . .   Subsequent  collisions  should  not  alter  this 
p      e: 

assertion. 

B.   FUNCTIONAL  DEPENDENCE  ON  OPTICAL  THICKNESS,  T 

P  monotonically  increases  with  increasing  t,  as  can  be 
clearly  observed  in  equation  (IV. 1).   The  rate  of  increase 
is  given  by  (using  equation  (IV. 1)) 


r~    v  v   ~  /      v  v? 

1-2  (-H)  sin  y  cos  6  +  (-S)  exp{n — H-r-  a  / 1-2  (-E)  sinycos6+  (-£)  } 
v      '       v     r  cos  y  V    v  v 


3P        , 
c       1 

3t     I  cos  y , 

m  m      '     '  v     m  m 


Here  the  effect  of  the  other  parameters  on  the  rate  with 
which  P   exponentially  approaches  one  is  explicitly  displayed. 
Appendix  E  contains  examples  of  the  collision  probability's 
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dependence  on  t  given  selected  values  for  the  speed  ratio, 
orbital  '  inclination ,  and  orbital  azimuth. 

The  physical  interpretation  of  this  functional  behavior 
is  simple.   Recall  that  for  measurements  at  appropriate 
wavelengths , 

x   =   /  Tir  n(r)hdr  . 

If  h  is  held  constant,  then  higher  values  of  t  imply  greater 
number  density,  n(r).   Intuitively,  one  would  expect  the 
collision  probability  to  increase  with  increasing  number 
density  (all  other  things  held  constant) ,  which  is  substan- 
tiated by  the  mathematical  expression  for  P  .   Similarly, 
if  n(r)  is  held  constant,  increasing  t  implies  that  h  must 
increase,  the  length  of  a  molecule's  path  must  increase, 
thus  increasing  its  exposure  to  collision. 

v 

C.   FUNCTIONAL  DEPENDENCE  ON  THE  SPEED  RATIO,  -£ 

v 

m 

The  behavior  of  the  collision  probability  as  a  function 
of  speed  ratio  is  very  interesting.   Figure  11  gives  a  repre- 
sentative example  of  the  variation  of  P   as  a  function  of 

v  /v   (see  Appendix  E  for  the  results  of  more  example  calcu- 
p  m 

lations.   For  6  less  than  yfT ,  a  minimum  for  P   occurs  at 

v  /v   greater  than  zero.   Why  does  this  happen? 

Some  insight  into  this  phenomenon  can  be  gained  by 

considering  the  partial  derivative  of  P   with  respect  to 

v  /v  ,  i.e., 
p  m 
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Figure  11.   Collision  Probability  as  a  Function  of  the 
Speed  Ratio  for  t  =  1.0  and  y  =  1/4tt  and 
Various  Values  of  9 
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With  all  variables  except  v  /v   held  constant,  P   has  a 

c        p'  m  c 

minimum  at 


v 

—2-  =   sin  y  cos  9  .  (IV. 3 

m 


From  Appendix  D,  the  right  hand  side  of  equation  (IV. 3)  is 
identically  equal  to  cos  a,  so  expressing  that  in  terms  of 
the  direction  cosine  gives 


v   =   v  cos  a  ;  (IV. 4) 


a  minimum  of  P   occurs  when  the  component  of  v   in  the  direc- 
c  m 

tion  of  v   is  equal  in  magnitude  to   v   .   Notice  that  relative 

p     *         *  '  p1 

to  variation  in  v  /v  ,  P   cannot  assume  a  minimum  value  for 

p  m   c 

speed  ratio  values  greater  than  one. 

Since  the  result  of  the  above  mathematical  analysis  may  . 
be  difficult  to  reconcile  with  intuition,  an  example  is 
given  to  clarify  its  physical  meaning.   Suppose  that  a  blind 
mouse  lives  in  a  hole  in  the  middle  of  a  single-lane,  one- 
way road.   Randomly  spaced  steamrollers  traverse  this  road, 
all  moving  with  speed  v  .   The  mouse,  smelling  cheese  in 
another  hole  a  few  meters  down  the  road,  dashes  for  that 

hole  with  speed  v   (in  the  same  direction  as  the  steamrollers). 

L  m 

Given  v  ,  for  what  v   will  the  mouse  have  the  best  chance 
m  p 

to  reach  the  cheese  unharmed,  neither  being  overtaken  by  a 
steam  roller,  nor  running  into  one?   The  answer  is  that  the 
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probability  of  collision  is  least  (in  this  case  zero)  when 

v   equals  v  .   For  v   greater  than  v  ,  steamrollers  will 
p   ^      m        p  ^  m 

tend  to  overtake  the  mouse;  as  v   gets  very  large  the  proba- 
bility of  collision  goes  to  one.   For  v   less  than  v  ,  the 

2  r  p  m' 

mouse  will  tend  to  run  into  a  steamroller.   The  mouse's 
situation  is  analogous  to  that  of  a  molecule  passing  through 
a  ring.   The  probability  of  collision  for  the  molecule  is 
minimized  when  it  "goes  with  the  flow"  of  the  particles  as 

much  as  is  possible.   This  minimum  value  of  P   will  not  in 

c  c 

general  be  equal  to  zero  for  molecules  passing  through  the 
ring,  as  it  is  for  the  mouse.   This  is  the  case,  because 
the  molecule  crosses  the  ring  plane  of  some  angle  (speci- 
fied by  y )  . 

A  final  comment  on  the  influence  of  v  /v   on  P   is  that 

p   m     c 


P 
c 


as 


v 
vm 


regardless  of  the  values  of  the  other  variables.   For  v 

P 

very  large,  the  ring  loses  its  discrete  character  becoming 
essentially  a  continuum j  collision  becomes  inevitable. 

D.   FUNCTIONAL  DEPENDENCE  ON  ORBITAL  INCLINATION,  y,    AND 
ORBITAL  AZIMUTH,  9 

y  and  6  specify  the  path  of  a  molecule  through  the 

ring.   Figures  12  and  13  show  two  examples  of  the 
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Fiqure  13.   P   as  a  Function  of  v  and  9  for  t  =  1 ,  — ^  =    2 
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functional  dependence  of  P   on  orbital  inclination  and 
azimuth. 

First,  consider  the  behavior  of  P  with  respect  to 

The  partial  derivative  with  respect  to  this  variable  is 


3Pc      t    v  Vd 

"5T  =  I  cos  y  I  ^T  sin^sine  ^-2  &  sinycose 
1    M  m  m 


v  ~  -1/2  /    v  v  - 

+  (-2)H    exp{n — =^-t- \  / 1-2  (-£)  sinycos9+ (-£)  }  . 
v        ^  cos  Y  V    v  v 


(IV.  5) 


From  equation  (IV. 5)  it  can  be  deduced  that  P   is  minimized 
when  -6  equals  zero  and  maximized  when  9  equals  tt   (when  all 
other  variables  are  held  constant).   As  expected,  both  Figures 
12  and  13  clearly  confirm  this  result.   Also,  notice  how 
rapidly  P   approaches  its  maximum  value  with  increasing  6. 
The  story  of  the  blind  mouse  from  Section  IV. C  provides  the 
physical  interpretation  for  this  behavior:   all  other  param- 
eters held  constant,  collision  probability  is  least  when 
going  with  the  flow  of  particles  (or  steamrollers) ,  and  is 
greatest  when  moving  against  them. 

The  variation  of  P   with  y  is  a  little  more  complicated. 
The  partial  derivative  of  P   with  respect  to  y  is 

8Pc  v  v     2  1/2 

-s —    =    [t  secy  tany  (1-2  (-*-)sin  ycos  6  +  (— )    ) 
3y  v  v 

m  m 


v  v  v     ~  -1/2  /         v  v     o 

-  x  (-£)cos9  (1-2  (-2)  sinycose-K-^)    )  ]  expin ^—r-  \  /1-2  (-£)  sinycosB-K-^)    } 

v  v  '  v  v    cos  y     V  v  '  v 

m  m  m  '  '  '    V  m  m 

(IV. 6) 
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From  equation  (IV. 6),  the  condition  for  minimum  P  ,  when 
all  variables  except  y  are  held  constant,  is 


v 
(-£)cos  6 

sec  Y ta^  Y    =    -^s_ —    .  (IV.7: 

1  +  2tan  y       L  +  (_p_}  2 

v 
m 

This  imposing  expression  has  a  simple  physical  meaning: 

(1)  if  v  /vm  is  less  than  one,  P   is  minimized  when 

the^magnitude  of  the  component  of  vm  in  the  direction 

of  v   is  equal  to  v  . 
P      M         P 

(2)  if  Vp/vm  is  greater  than  one,  Pc  is  minimized  when 
the^magnitude  of  the  component  of  vp  in  the  direction 
of  vm  is  equal  to  vm. 

Once  again,  the  blind  mouse  analogy  may  be  helpful  in  under- 
standing this  behavior. 

Before  leaving  this  section,  a  few  words  are  in  order 
to  clarify  possible  concern  over  the  point  of  "undefined" 
collision  probability  shown  in  Figure  12.   This  anomaly  only 

occurs  when  v  /v   is  equal  to  one.   For  this  particular  speed 

p   m     ^  r  ^ 

ratio,  when  9  is  zero,  and  y    is  tt/2,  then  the  molecule  is 
moving  in  concert  with  the  ring  particles.   Thus,  if  it  is 
not  initially  in  contact  with  a  particle,  it  never  will  be. 
On  the  other  hand,  if  initially  in  contact  with  a  particle, 
it  will  never  separate  from  that  particle.   Therefore,  the 
notion  of  a  collision  probability  is  meaningless  under  this 
unique  circumstance. 
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V.   DISCUSSION 

A.  GENERAL  COMMENTS 

To  effectively  exploit  the  collision  theory  developed  in 
this  paper,  one  generally  must  couple  it  with  a  ballistic 
transport  model  which  provides  values  for  the  speed  ratio, 
orbital  inclination,  and  orbital  azimuth.   More  accurate 
estimation  of  atmospheric  density  will  result,  although 
improvements  in  this  area  are  limited  by  the  lack  of  knowledge 
of  details  of  a  molecule-particle  collision  (i.e.,  the  stick- 
ing coefficient  described  in  Section  I) .   The  frequency  of 
collisions  also  affects  the  spatial  extent  of  the  ring 
atmosphere.   The  trajectory  of  a  molecule  is  initially 
determined  by  the  nature  of  the  production  process,  while  its 
evolution  is  dictated  by  collisions. 

Consideration  of  the  special  case  of  low  energy,  iso- 
tropic production  processes  (e.g.,  sublimation)  allows  one 
to  apply  the  collision  theory  in  making  inferences  on  the 
spatial  extent  of  the  ring  atmosphere. 

B.  APPLICATION  TO  AN  ISOTROPIC,  LOW  ENERGY  PRODUCTION 
PROCESS 

Define  an  isotropic,  low  energy  production  process  to 

be  a  mechanism  which  ejects  molecules  from  ring  particles 

such  that: 

(1)   the  direction  of  ejection  relative  to  the  particle 
is  random. 
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(2)   v    <<  v  ,  which  implies  v  /v   =  1. 
ej  p  r      p   m 

The  orbit  of  a  molecule,  prior  to  collisions  with  ring 
particles,  will  be  very  nearly  circular  with  a  radius  equal 
to  that  of  its  parent  ring  particle.   The  orbital  plane 
will  be  only  slightly  inclined  with  respect  to  the  ring 
plane,  e.g.,  for  sublimation  the  orbital  plane  is  inclined 
at  most  1°  (see  Section  IV. A) . 

When  a  collision  does  occur,  two  outcomes  are  possible: 
the  molecule  "sticks,"  thus  removing  it  from  the  atmosphere; 
or  the  molecule  rebounds,  assuming  a  new  trajectory.   In 
the  latter  case,  what  can  be  said  about  the  new  orbit? 
Since.,  the  collision  will  surely  be  inelastic  to  some  degree, 
the  speed  of  the  rebounding  particle  relative  to  the  ring 
particle  is  less  than  its  incident  speed.   The  speed  ratio 
remains  almost  equal  to  one,  implying  a  circular  orbit. 
Following  the  approach  used  in  Section  IV. A,  especially 
Figure  10,  the  primary  effect  of  the  inelastic  collision  is 
to  reduce  the  inclination  of  the  orbital  plane.   Thus, 
collisions  act  to  reduce  the  spatial  extent  of  the  ring 
atmosphere;  the  details  of  this  reduction  being  determined 
by  the  frequency  of  collisions,  and  the  actual  inelasticity 
of  an  individual  collision. 

The  more  frequent  collisions  are,  the  more  closely  con- 
fined to  the  vicinity  of  the  ring  will  be  the  ring  atmosphere 
The  collision  theory  can  be  used  to  examine  collision  fre- 
quency.  Specifically,  the  average  probability,  as  a  function 
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of  time,  that  a  molecule  has  not  suffered  at  least  one 
collision  is  calculated  below. 

When  a  molecule  is  initially  ejected,  it  may  collide 
with  a  ring  particle  before  ever  leaving  the  ring  plane. 
To  calculate  the  probability  that  it  collides  with  at  least 
one  ring  particle,  adopt  a  frame  of  reference  in  which 
the  particles  are  fixed.   The  probability  that  the  molecule 
is  ejected  from  a  particle  in  dz '  is  (see  Figure  14) 


P(z' )dz'   = 


h   dz' 
h 


0  <  z '  <  h 


(V.l) 


while  the  probability  that  it  is  ejected  into  solid  angle 
dQ.    is 


P(ft)dft   = 


sin  y  '  dy  '  ,   0  <  y '  <    ^    i 


(V.2) 


"Effective 
Thickness" 


£. 


Ejection  Angle 


!      Ejection  Location 


Figure  14.   Ejection  of  a  Molecule  from  a  Particle  at  z1 
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since  the  production  process  is  isotropic.  The  probability 
density  for  a  molecule  ejected  at  y'  from  z*  colliding  with 
a  ring  particle  before  leaving  the  ring 


,       r  -t  (h-z  ' )  !  „      .    Tf 

1  -expir-— - r1-}  i       0  <  Y   <  o- 

^  h  cosy '  '     —  '     2 

Pc( 


i         r   +T-Z  '   i        7T       ,. 

r  h  cos  y  '       2  —  '   — 


which  follows  directly  from  equation  (III. 12),  when  the 
effective  slab  thickness  is  taken  into  account.   The  aver- 
age probability  of  at  least  one  collision  for  a  molecule 
produced  by  an  isotropic  process  is  found  by  integrating  over 
the  product  of  equations  (V.l),  (V.2),  and  (V.3): 


2"1   h 


Pc  =  J   /  |hsin  V  u-«p{=|JfcS^.]d«'dy 


+  J*  lo   |hsin  Y.[l-exp{jr^jr.}]dz'  dY'  . 

2n 


Making  a  change  of  variables,  x  =  cos  y'  and  x  =  -  cosy1 
in  the  first  and  second  integrals,  respectively,  then  inte- 
grating, yields 


P    =   1  -  -  (      x(l  -  e  T/x)dx  (V.4 

C  T  J  n 
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Therefore,  the  average  probability  that  a  molecule  initially 
(i.e.  ,  t  =  0)  escapes  the  ring  without  collision  is 


esc 


=   1  -  P 


1  f         ,,     -t/xx  . 
-J      x  (1  -  e     )dx 


:v.5 


Table  2  gives  P     for  values  of  t  of  interest  for  the  rings. 
r       ecs 

Notice  that  there  is  a  significant  chance  of  collision  immedi- 
ately after  a  molecule  is  ejected;  suggesting  that  the  spatial 
extent  of  the  atmosphere  is  biased  toward  being  closely 
confined  to  the  ring  even  before  subsequent  evolution  of  the 
system  is  considered. 

TABLE  2 

AVERAGE  PROBABILITY  THAT  A  MOLECULE  INITIALLY 

(t  =  0)  ESCAPES  THE  RING  PLANE  WITHOUT 

COLLIDING  WITH  RING  PARTICLES 


T 

esc 

0.1 

0.837 

0.5 

0.556 

1.0 

0.390 

1.5 

0.296 

2.0 

0.235 

During  each  orbital  period,  a  molecule  in  a  circular 
orbit  traverses  the  ring  two  times.   If  that  molecule  has 
not  collided  with  a  ring  particle,  y1  retains  its  initial 
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value  for  each  subsequent  passage  of  the  rings,  and  the 

probability  that  it  does  not  collide  with  a  ring  particle 

-2t   .  14 

during  each  period  is  exptn ^—tts  .     The  probability  that 

|  cos    y     I  r  -1 

the   molecule    is    ejected    into    solid   angle   dft  without   colliding 
with   a   ring   particle    is: 


1  h  1  1 

jsiny't/     jrd-P    (y' ,z')  )dz']dy'     =  ^siny1. 


E2£Ll(i-exp{— '-?—}),   0  <y'  <tr 
t  *  cosy'      '       — '        2 


-COSY     /i  r      T      -i  x       1      ^     .    ^ 

— (1-expl r}    ,  ^-tt  <y'  <r> 

x  c  cosy  2    —      — 


If  T-  is  the  number  of  completed  orbital  periods,  then  the 
p 

probability  that  the  molecule  has  not  collided  with  any 
particle  as  a  function  of  time,  given  in  terms  of  T  ,  is 


C2£Y_!_(i_exp{— 'ST}t   o  <Y'  <iir 
t      ^  cosy     -    2 


TT 


-2tT 


P   (T  )  =    isiny1/  >exp(-, — .-rdy'  , 

esc  p     '  n  2   '  \  /  v   |cos  y  I 


0 


-cosy1  ,,     r  t  i  1    , 

— T-(1-exP{E3i7r}'2TT^Y  ^ 


or 


,   1    -2tT  /x       7 
P    (T  )   =   -  /   xe     p   (1-e  T/X)dx  .  (V.6) 

esc   p      T   0 


14 

This  result  can  be  easily  understood  by  considering 

a  slab  of  thickness  2h.   Extending  this  reasoning,  if  Tp 

periods  have  been  completed,  the  probability  of  no  collision 

for  a  molecule  initially  escaping  is  exp{ -2tT  / | cosy ' | } . 


62 


Table  3  gives  the  probability  that  a  molecule  has  remained 
free  of  collisions , as  a  function  of  time  (given  in  orbital 
periods)  for  various  optical  thicknesses.   These  results 

TABLE  3 

VARIATION  OF  P     WITH  T  AND  T 
esc  p 


esc 


T 

P 

T    =    0.1 

t    =    0.5 

T    =    1.0 

T    =    1.5 

t    =    2.0 

0 

0.837 

0.556 

0.390 

0.296 

0.235 

1 

0.519 

0.106 

0.021 

0.005 

0.001 

2 

0.357 

0.028 

0.002 

1.56xl0"4 

-o 

3 

0.255 

0.008 

1.94xl0"4 

~0 

-o 

4 

0.186 

0.002 

-0 

-0 

-o 

5 

0.138 

0.001 

~0 

~o 

~o 

strongly  imply  collisions  are  quite  frequent  for  isotropic, 
low  energy  production  processes.   This  in  turn  implies  that 
the  ring  atmosphere  associated  with  this  process  will  tend 
to  be  more  concentrated  in  the  vicinity  of  the  ring  than 
analysis  of  the  distribution  of  particles  prior  to  collision 
would  indicate. 

C.   ON  MULTIPLE  COLLISIONS 

The  collision  probability  theory  developed  in  this  thesis 
cannot  be  used  to  precisely  calculate  the  probability  of 
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multiple  collisions,  i.e.,  the  chance  that  a  molecule  col- 
lides with  two,  three,  etc.,  particles  during  a  single 
passage  through  the  rings.   Computer-simulation  models  may 
provide  the  best  way  to  approach  this  problem. 

Based  on  the  results  of  the  preceding  section,  the 
probability  of  at  least  one  collision  is  quite  high  in  the 
A  and  B  rings.   This  suggests  that  multiple  collisions  could 
be  a  significant  factor  in  analysis  of  the  ring  atmosphere. 
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VI.   CONCLUSIONS 

In  this  paper,  a  model  for  collisions  between  ring 
particles  and  molecules  in  the  ring  atmosphere  has  been 
developed.   The  resulting  expressions  for  the  probability 
of  at  least  one  collision  when  a  molecule  passes  through 
the  ring  are 


P    =   1  -  exp{n "T  ,i  }  ,  (III.  12 

c  r  I  cos  y  I 


for  a  frame  of  reference  where  the  ring  particles  are  fixed, 
and 


Pc  =  l-exp{|cQ~T  ,  \/l-2(jft  sin y  cose +(^)2}  ,  (IV.l) 

1    i  I  v     m  m 


for  a  frame  of  reference  fixed  on  Saturn  in  which  the  parti- 
cles are  seen  to  orbit  the  planet.   They  are  more  complete 
than  those  currently  employed  in  analysis  of  the  ring  atmos- 
phere, since  specifics  on  molecular  trajectory  and  the  rela- 
tive velocity  of  a  molecule  with  respect  to  the  ring  parti- 
cles are  included.   Thus,  this  analytic  collision  theory  is 
compatible  with  computer-simulated  ballistic  transport 
models,  since  they  include  similar  effects. 

The  frequency  of  collisions  influences  both  the  total 
population  of  molecules  in  the  ring  atmosphere,  and  their 
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spatial  distribution,  as  shown  by  Dennefeld  (1974).   Since 
knowledge  of  molecular  sources  is  still  evolving,  and  under- 
standing of  the  dynamics  of  an  individual  collision  (e.g., 
Dennefeld ' s  "sticking  coefficient")  is  in  its  infancy, 
improved  estimates  of  molecular  populations  cannot  be  made 
based  on  this  revision  of  collision  theory,  alone.   Infer- 
ences on  spatial  distribution  can,  however,  be  made.   When 
a  low  energy,  isotropic  production  process  is  considered  for 
optical  thicknesses  characteristic  of  the  A  and  B  rings  of 
Saturn,  collisions  are  found  to  be  quite  frequent.   This 
implies  an  atmosphere  closely  confined  to  the  vicinity  of  the 
rings..   Extending  this  result  to  more  energetic  production 
processes  (e.g.,  see  Section  III. A)  implies  that  the  ring 
atmosphere  will  be  toroidal  in  shape.   Thus,  assuming  a 
toroidal  shape  will  probably  be  more  productive  than  assum- 
ing spherical  shape  in  future  studies  of  the  ring  atmosphere. 

As  a  final  note,  the  collision  theory  used  here  in  con- 
nection with  the  ring  atmosphere  should  be  useful  in  another 
related  application.   Ballistic  transport  of  macroscopic 
particles  produced  during  meteroid  impacts  may  account  for 
some  of  the  radial  structure  in  Saturn's  rings  (Ip,  19  83; 
Durisen,  1984).   These  particles  form  a  "chip  halo"  in  the 
vicinity  of  the  rings.   As  long  as  the  chips  are  small  rela- 
tive to  ring  particles,  the  model  developed  here  applies 
for  calculation  of  collision  probabilities. 
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APPENDIX  A 
RING  NOMENCLATURE  AND  DIMENSIONS 

Saturn's  ring  system  consists  of  seven  distinct  rings. 
The  A,  B,  C,  D,  and  F  rings  are  characterized  by  distinct 
boundaries.   Their  nominal  thicknesses  are  not  more  than 
-200  m  (Cuzzi  et  al.,  1984),  and  may  well  be  less  than  5-10 
m  (Bridges  et  al . ,  1984) .   The  E  and  G  rings  are  diffuse 

structures  with  poorly  defined  boundaries  and  nominal  thick- 

3        2 
nesses  of  ~10   and  10   km,  respectively  (Cuzzi  et  al . , 

19  84^.   Table  Al  is  a  summary  of  the  radial  structure  of  the 

ring  system.   Note  that  the  F  and  G  rings  are  quite  narrow. 

TABLE  Al 
RADIAL  STRUCTURE  OF  SATURN'S  RINGS 

Boundaries 

Ring/Region         R                  km  mass 

D           -1.11-1.235      66,970-74,510  ? 

C            1.235-1.525     74,510-92,000  2  x 10_9  Ms 

B            1.525-1.949     9  2,000-117,58  0  5  x 10_8  Ms 

Cassini         1.949-2.025     117,580-122,170     1  x 10~9  Ms 
Division  '        ' 

A  2.025-2.267     122,170-136,780     1  x lo"8  Ms 

F  2.324  140,180  ? 

G  2.82  170,100  1   xio"17    Ms 

E  3-8  181,000-438,000  ? 

2  6 
R   =  60,330  km  (Saturn's  radius)   M   =  5.685x10    kg  (Saturn's 
s  s  \ 

mass) 
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The  Cassini  Division  is  a  well-defined  gap  between  the  A  and 
B  rings,  which  is  by  no  means  "empty." 

The  optical  depth,  t,  is  a  dimensionless  parameter 
directly  related  to  mass  density.   Table  A2  gives  average 
optical  thickness  for  the  main  rings.   Optical  thickness  in 
the  other  rings  is  much  lower.   Since  optical  thickness 
increases  with  increasing  mass  density,  it  follows  that  the 
A  and  B  rings  should  be  the  most  significant  sources  of  a 
ring  atmosphere. 

TABLE  A 2 

.AVERAGE  OPTICAL  THICKNESS  OF  THE  MAIN  RINGS  BASED 

ON  STELLAR  OCCULTATION  OBSERVATIONS  FROM 

VOYAGER  2  (ESPOSITO  ET  AL . ,  1984) 

REGION              BOUNDARIES  (R  )              TAVE 
s_  

inner  C  1.24-1.39  0.08. 

outer  C  1.39-1.52  0.15 

inner  B  1.52-1.66  1.21 

middle  B  1.66-1.72  1.76 

outer  B  1.72-1.95  1.84 

Cassini  Division  1.95-2.02  0.12 

inner  A  2.02-2.16  0.70 

outer  A  2.16-2.27  0.57 


R    =   60,330  km 
s 
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APPENDIX  B 

DERIVATION  OF  P   WHEN  v   IS  MOT  EQUAL  TO  ZERO 
c p a 

In  Section  II  of  the  main  body  of  this  paper,  an  expres- 
sion for  the  probability  that  a  molecule  passing  through  a 
ring  collides  with  at  least  one  ring  particle,  P  ,  was 
derived  in  a  frame  of  reference  in  which  the  ring  particles 

are  fixed  (i.e.,  v   equal  to  zero) .   Here  it  will  be  demon- 

P 

strated  that  the  same  expression  for  P   results  when  the 
derivation  is  made  in  a  frame  of  reference  where  the  slab 
translates  with  v   not  equal  to  zero. 

The  model  upon  which  this  derivation  is  based  is  the 
same  as  that  described  in  Section  II.   Reviewing,  a  molecule 
passes  through  a  region  in  space  which  contains  randomly 
distributed  (spatially)  spherical  particles.   This  region, 
the  "slab,"  is  of  large  extent,  L,  in  two  dimensions-  and  of 
relatively  small  thickness,  h.   The  entire  slab  moves  through 

space  with  constant  velocity,  v  ,  in  a  direction  normal  to 

->• 
the  slab's  thickness  (v   is  measured  relative  to  some  iner- 

P 

tial  frame  of  reference).   The  radius,  or  "size,"  of  each 

particle  is  restricted  to  an  interval  [r    ,r    ]  with  r 

c  min   max         mm 

much  greater  than  the  nominal  radius  of  the  molecule,  and 

with  r    much  smaller  than  h.   The  size  distribution  of  the 
max 

particles  may  be  either  discrete  or  continuous.   The  number 
of  particles  in  the  slab  is  further  assumed  to  be  large. 
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Finally,  the  number  density  as  a  function  of  size  is  denoted 
by  n(r),  defined  as  the  number  of  particles: 

(1)  of  radius  r  per  unit  volume  (discrete  size 
distributions) . 

(2)  in  the  radius  interval  [r,r+dr]  per  unit  volume 
(continuous  size  distributions) . 

What  is  the  probability  that  the  molecule  collides  with 
a  particle  in  the  slab?   In  Section  II  this  question  is 
answered  with  relative  ease  by  converting  to  a  frame  of 
reference  in  which  the  particles  are  at  rest.   Below,  no 
such  simplification  is  made. 

When  v   is  not  equal  to  zero,  the  notion  of  cylindrical 
"collision  regions,"  used  in  Section  II,  loses  its  validity. 
Another  approach  must  be  taken.   Suppose  that  the  slab  is 

subdivided  into  regions  of  very  small  thickness,  Ah,  and 

2 
finite  area  L   (see  Figure  Bl).   If  the  probability  that  the 


/\ 


h 


X  Ah. 


f\ 


Figure  Bl .   Edgeview  of  the  Slab  and  a  Representative 
Subregion 
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molecule  does  not  collide  with  particles  of  size  r  whose 

centers  lie  in  a  specific  sublayer  can  be  calculated  (call 

this  P   (r,Ah.)),  then  the  total  probability  of  no  colli- 
no      1  c  2 

sion  with  particles  of  size  r  will  be 


Pno(r>   "   n  Pno(r'4hi»  '  (B-1) 

1 

since  all  of  the  small  sublayers  are  mutually  independent. 
The  collision  probability,  following  Section  II,  is 


P    =   1  -n  P   (r.)   =   1-nriP   (r . ,Ah. )  (B.2) 

c  no  -\  no  -]        i 


for  discretely  distributed  r. 

The  problem  now  is  to  find  an  expression  for  P   (r.,Ah.) . 
1  no   3    i 

First,  the  case  of  discrete  r  is  solved  from  which  the  con- 
tinuous solution  can  be  deduced.   Adopt  a  point  of  view  in 
which  the  molecule  has  radius  r,  and  the  particles  are  point 
masses  (on  reflection  this  is  analogous  to  the  real  situation 
as  far  as  collision  probabilities  are  concerned).   Then, 

P   (r.,Ah.)  can  clearly  be  seen  to  be  equivalent  to  the 
no   i    i  2  ^ 

probability  that  the  volume  in  the  subslab  "swept  out"  by 


Two  observations  on  equation  (B.l)  are  in  order.   First, 
there  are  infinitely  many  factors  in  the  product  over  i. 
Second,  Pno(r,Ahj_)  must  be  the  same  for  all  i,  since  in  this 
model  the  statistical  properties  of  each  subregion  are  iden- 
tical.  Equation  (B.l)  may  be  difficult  to  interpret  at 
this  point,  but  its  meaning  will  be  clarified  later  in  this 
appendix. 
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the  passing  molecule  has  no  particles  centered  within  it 
The  probability  of  no  collisions  with  particles  of  size 
r.  in  subslab  Ah.  is: 


L  Ah. -(Vol.  Swept  Out)  n(r.)L  Ah. 
Pno(r  ,ah)  =  [ i--- ]   3     1     ,  (B.3 

J  L  Ah. 


since  the  particles  are  assumed  to  be  randomly  distributed 
throughout  the  slab  (see  Section  II  for  a  more  detailed 
explanation  of  this  assertion).   In  general,  the  mathematics 
involved  in  finding  this  "volume  swept  out"  is  extremely 
tedious.   To  demonstrate  the  method,  the  special  case  of  y 
equal  to  zero  or  tt  (molecular  path  normal  to  the  plane  of 
the  slab)  is  solved  below. 

On  what  portion  of  its  path  can  a  molecule  potentially 
collide  with  a  ring  particle  of  size  r.  in  Ah . ?   Consideration 
of  Figure  B2  shows  that  the  pertinent  part  of  the  path  is  of 
length  2r . .   If  time,  t,  is  set  equal  to  zero  when  the  mole- 
cule initially  encounters  the  subslab  (Figure  B2(a)),  then 
at  any  t  such  that 


2r  . 

0   <   t   < 


—   v 

m 


the  molecule  will  collide  with  any  particle  is  a  disk  of 
radius  (call  this  the  "interaction  radius") 

2  2  ^2 
RT(t)   =   (2rv  t  -  v  t  )     ,  (B.4) 

I  mm 
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Molecular 
Path 


Ah. 

1 


T 


t  >  0 
(b) 


Note : 
Y  =  0 


Figure  B2. 


Relationship  Between  the  Molecule  and  a 

Subslab  When  Collision  with  Ring  Particles 

Centered  in  the  Subslab  is:   (a)  First 

Possible  (t  =  0) ;  (b)  At  Any  Time,  t, 

Such  That  0  <  t  <  2r./v 

J     m 


where  R   is.  defined  as  shown  in  Figure  B2(b)  .   As  the  inter- 
action radius  evolves,  the  particles  in  the  subslab  move 
with  speed  v  .   Considering  these  two  factors  together,  it 
can  be  deduced  that  an  oblong  "interaction  region"  in  the 
subslab  must  be  particle-free,  if  no  collisions  are  to  occur. 
To  actually  calculate  this  volume,  first  relate  R  to  the 
distance  that  the  slab  moves  in  time  t,  call  this  D  ;  clearly 


D_(t)   =   v  t  . 

I  P 


B.5 


Substituting  from  equation  (B.5)  into  equation  (B.4)  yields 


v        v   ~    1/2 
RI(DI)   =   (2r(^L)DI  -  (^Vd^) 


(B.6 


73 


Figure  B3  is  a  sketch  of  equation  (B.6) .   The  oblong  area 

in  Figure  B3   (bounded  by  the  graph  for  R   and  its  reflection 

across  the  D  -axis,  which  is  shown  by  a  dashed  line)  cannot 


RI(DI) 


*  D. 


Figure  B3.   R   as  a  Function  of  D 


represent  the  projection  of  the  interaction  region  into  the 
plane  of  the  subslab.   Consideration  of  instantaneous  inter- 
action region,  represented  by  the  circle  in  Figure  B3 ,  illus- 
trates that  the  actual  interaction  region  is  larger.   Knowing 
R  (D  ),  the  function  enveloping  the  interaction  region  can 
be  found.   Call  this  function  R  . 

Consider  Figure  B4 ;  it  shows  one-quarter  of  the  inter- 
action volume's  projection  into  the  plane  of  Ah. .   At  a  given 

time  t",  D'  equals  v  t1,  and  the  interaction  radius  is  R ' . 
I   ^       m  I 

The  value  of  R" (D  )  for  any  D   is 


Rj(Dj.) 


2         2  1/21 
([R^Dj)  [  -[D'-DI]Z)1/ZJ 


Dj-DJI  <R^(D^) 


otherwise  . 


(B.7) 
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VRi 


Figure  B4 .   Geometry  for  Calculating  r! 


D. 


R  (D  )  is  the  maximum  value  of  R"(D|).   Setting  the  first 
derivative  of  R"  with  respect  to  Dl  equal  to  zero,  and 
solving  for  Dj  gives 


4 


J   P 

[1  +  (^)2] 
P 


(B.8) 


when  equation  (B.6)  is  substituted  into  equation  (B.7). 
Substituting  from  equation  (B.8)  into  equation  (B.7)  gives 
the  maximum  value  of  R"(D|)  for  a  given  D  ,  i.e., 


VV 


[c£h 

P 


n  V  -  V 

Z  .    m>  ,  ~   ^    ^2  ,  m> 
r  (   ,  +2rjD  -DC—)   1/2 

_i_£ { p_}] 

V   «, 

1  +  (^)2 

V 

p 


(B.9) 
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Using  equation  (B.9),  the  region  in  the  sublayer  that  must 
be  devoid  of  particles  of  size  r.  if  no  collision  is  to 
occur  between  the  molecule  and  such  particles,  has 

v 

r.(-P) 

3  vm 

Volume   =   (Ah.)  4/        R*(D  )  dD   ,  (B.10) 

1     D  ^ 

I (min) 


where 


DI(min)      "       [£>-U+    <^>2>1/2)r      •  (B.ll) 

m  m        J 


When  the  integral  in  equation  (B.10)  is  evaluated, 


7  Vr,     1      X/2 

Volume   =   irr  {1  +  (-£•)}    Ah .  ,  (B.12 

J  m 


and  substituting  into  equation  (B.3)  gives 


?       o     v   ?  1/2  7 

LzAh.-7Tr f{l+(-£)  }    Ah.  n(r.)L  Ah. 
l    j     v  1     J      1 

P   (r.,Ah.)   =  [ 5 5} ]  •    (B.13a! 


Equation  (B.13)  can  be  rewritten  as 

7rrZn(r.)  {1  +  (-£)  }     „,,.  ,T2.. 
j    j       v         n  (r  .  )  L  Ah . 

P   (r.,Ah.)   =  [1 5 2 ]    ^      \      (B.13b 

n°  1         x  L2n(r.) 
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Summing  over  i  yields  the  probability  that  no  collision 
occurs  in  the  entire  salb,  i.e., 


9  v   -  1/2 

7TrZn(r  .)  {1  +  (-£)  }      ,   .  T  2. 
:    J      v        n(r.)L  h 

P„„(r.)   =   [1 5 - ]    ]  (B.14a 

n°   J  L2n(r.) 


Multiplying  the  second  term  in  the  brackets  in  equation 
(B.14)-  by  h/h  yields 


v   ,  1/2 

,r2n(r.){1+(_£)  }    h  n(r  )L2h 

Pno(r  )   =   [1 m ]    3  (B.14b 

J  L  n  (r  .  )  h 

2 
but  the  total  number  of  particles  of  size  r.,  n(r.)L  h,  is 

assumed  to  be  large,  so  equation  (B.14b)  can  be  rewritten 

in  exponential  form  as 


v   ~  1/2 

Pno(r.)   =   exp{-Tir  n(r  .  )  {1  +  (-H)  }    h}  (B.14c) 

-1  -1    ■*  vm 


Therefore,  the  collision  probability  is  (from  equation  (B.2 


7  v  7    1/2 

P    =   1  -  n  exp{-Trr  n(r  .)  {1  +  (-£)  }    h}  , 
c  1    1       v 

j         J    J        m 


P  v   ?  1/2 

P^   =   1  -  exp{-rr(J  rTn(r.)){l  +  (-£)  }    h}  ,  (B.15a) 

c  j   D    D       vm 


for  r.  discrete,  or  by  extension 
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7  Vn  ■)     l/2 

P    =   1  -  exp{-TT(/  r  n(r)dr)  {1  +  (-£)  }    h}  ,         (B.15b) 

w  V 

m 


for  continuous  r. 

Comparison  of  equations  (B.15a)  and  (B.15b)  with  equations 
(II. 7)  and  (II. 8)  (with  y    equal  to  zero),  reveals  the  equiva- 
lence of  solving  the  collision  probability  problem  in  two 
different  frames-of-ref erence.   That  the  results  are  the  same 
is  not  surprising.   What  is  interesting  is  the  relative 
difficulty  of  the  two  approaches. 
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APPENDIX  C 
FRAME  OF  REFERENCE  CONVERSION 

In  Section  II  of  this  thesis,  the  collision  probability, 

P  ,  is  calculated  in  a  frame  of  reference  in  which  the  ring 

c  ' 

particle  velocity,  v  ,  is  equal  to  zero.   When  coupling  the 
collision  theory  with  ballistic  transport  models,  a  more 
convenient  frame  of  reference  is  a  localized  cartesian  sys- 
tem fixed  in  the  vicinity  of  the  intersection  of  the  mole- 
cule's path  and  the  ring.   The  x-axis  is  in  the  direction  of 
motion  of  the  ring  particles;  the  y-axis  is  directed  radially 
outward  from  Saturn;  and  the  z-axis  is  normal  to  the  ring 
plance  (see  Figure  CI).   Here  the  direction  cosines  of  a,  (3 , 


Y 


x 


Figure  CI.   A  Coordinate  System  Fixed  at  a  Specific 
Location  in  the  Rings 
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and  y  quite  accurately  describe  the  path  of  the  molecule 
through  the  ring  which  is  essentially  linear  owing  to  the 
great  radius  of  curvature  of  the  molecular  trajectory 
(-10   km)  relative  to  the  extreme  thinness  (-10    km)  of  the 
rings . 

Let  the  frame  of  reference,  for  which  the  ring  particles 
in  the  region  through  which  the  molecule  passes  are  at  rest, 
be  called  the  "primed"  system.   The  simplest  coordinate 
system  that  can  be  used  to  describe  the  molecule's  path  in 
this  frame  of  reference  will  have  the  x1-,  y'-,  and  z'-axes 
in  the  same  direction  as  the  corresponding  axes  in  the  un- 
primed  system  with  the  y ' -axis  in  the  direction  of  v   (see 
Figure  C2).   a',  B',  and  y'  yield  direction  cosines  describing 
the  molecules  path  in  the  primed  system. 


x 


V 


Figure  C2.   Direction  Cosines  in  the  Two  Coordinate  Systems 


In  Section  II  the  collision  problem  is  attacked  using 
the  primed  system.   The  expressions  for  P   derived  there 
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(equations  (II.  5)  and  (II.  6))  contain  a  factor  of  j cos  y'j 

An  expression  for  this  factor  in  terms  of  a,  B,  y,  v  •,  and 

v  must  be  found  to  represent  P   in  the  unprimed  system, 
p  r        c  r       2 

First,  note  that: 


(C.l) 

m      m    p 


gives  the  relationship  between  the  molecular  velocities  in 
the  two  systems.   Writing  vector  equation  (C.l)  in  component 
form  yields: 


v1    =v-v;  ( C .  2  a 

mx      mx     p 


v1    =   v    ;  (C.2b 

my      my 


v1    =   v  (C.2c 

mz      mz 


keeping  in  mind  the  definitions  of  the  two  systems  given  in 
Figure  CI.   Equations  (C.2)  can  be  rewritten  in  terms  of 
the  direction  cosines: 


v'  cos  a'   =   v   cos  a  -  v   ;  (C.3a 

m  m  p 


v1  cos  S1   =   v   cos  8  ;  (C.3b 

m  m 


v'  cos  Y1   =   v   cos  Y  •  (C.3c 

m      '        m      ' 
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Squaring  equations  (C.3)  and  adding  yields 


22        2        2  2222 

v'  (cos  a'  +cos  3'  +cos  y')      =   v  (cos  a  +cos  S  +cos  y) 

2 

+   v   -  2v  v  cos  a  .       (C.4 
P     m  p 


Since  the  sum  of  the  squares  of  direction  cosines  is  iden- 
tically equal  to  one,  equation  (C.4)  becomes 


,2       2 

v   + 


v'    =   v   +v   -2vv   cos  a  .  (C.5 

m        m     p      m  p 


From  equation  (C.3c) 


cos  y 
m      'm  'cos  y 


v:    =    v^(_:  ;.)  ,  (c.6) 


which,  when  substituted  into  equation  (C.5),  yields 


2  ,cos  y  \  2       2,2    -  ,_  _ 

v  ( —)         =   v   +v   -2vv   cos  a  .         C.7 

m  cos  y         m    p      m  p 


Solving  for  |cos  y' I    gives 


,        /      V  V    9 

— rA/i-2(-£)cos  a  +  (-£)   .      (C.8) 

:os  y  V     v  v 

1  '  V      m  m 


cos  y ' , 

m         m 


Equation  (C.8)  is  the  fundamental  result  of  this  derivation. 
When  substituted  into  equations  (II. 5)  and  (II.  6),  the  result- 
ing expressions  for  P  ,  equations  (II. 7)  and  (II. 8),  have 
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the  advantage  that  they  are  written  in  terms  of  the  frame 
of  reference  most  often  utilized  when  studying  the  rings. 
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APPENDIX  D 
RELATING  ORBITAL  AZIMUTH  TO  THE  DIRECTION  COSINES 


Consider  a  unit  vector,  u,  as  shown  in  Figure  Dl .   Its 
direction  can  be  specified  by  the  direction  cosines  (a,B/Y)/ 
or  by  (y,9)  where  9  is  called  the  "orbital  azimuth."   How 
are  these  two  representations  related? 


>    Y 


Figure  Dl .   Angles  Used  to  Specify  the  Direction  of  u 


First,  note  that 


2        2        2 
cos  a  +  cos  6  +  cos  y      =      1  / 


or 


2fl      ,       2        2 
cos  3   =   1  -  cos  a  -  cos  y 


(D.l) 
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The  x  and  y  components  of  u  are 


u    =   cos  a    ,  (D.2a) 

2\. 


u   =   cos  6  .  (D.2b) 

y 


By  definition 


u 

tan  9   =   -^  ,  (D.3 

u 
x 


and  substituting  from  equations  (D.2)  gives 


tan  0   =   S2S-1  .  (D.4) 

cos  a 


Squaring  equation  (D.4)  and  substituting  from  equation  (D.l) 
yields 


2        2 
,   2_      1  -  cos  a  -  cos  Y  /-n.    c 

tan  6   =   2 •  (D.5 

cos  a 


Rearranging  equation  (D.5)  gives 


2      2  2 

(1  +  tan  9)cos  a   =   1  -  cos  y    , 


but 


2  2  2  2 

1  +  tan  9   =   sec  9     and     1  -  cos  y      =   sin  y    , 
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so 


2a     2  2 

sec  9  cos  a   =   sin  y 


(D.6) 


Finally, 


cos  a   =   sin  y  cos 


(D.7) 


which  is  useful  for  expressing  Pc  in  terms  of  y    and  0. 
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APPENDIX  E 

SOME  RESULTS  OF  CALCULATION  OF  P 
c 

This  appendix  consists  of  several  examples  of  results 
for  the  calculation  of  P  ,  based  on  the  equation 


v  v 


P    =   1  -exp{,    T  ,Wl-2(^)sin  ycos  9  +  (-£)  2  } 
^  I t-u&  Y  I  V      m  vm 


Calculations  were  made  over  the  range  of  optical  thickness 
pertinent  to  Saturn's  A  and  B  rings: 

0.1  <_     t   <   2.0  . 

The  figures  have  been  arranged  to  allow  one  to  compare  the 
effect  of  varying  the  parameters  of  the  equation. 
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Figure  El.   Collision  Probability  for  y    =  0 
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Figure   E2 .      Collision   Probability    for    y    =   =■ 
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Figure  E4 .   Collision  Probability  fo.r  y    =  p-  tt 
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APPENDIX  F 
ANNOTATED  BIBLIOGRAPHY 

The  purpose  of  this  appendix  is  to  briefly  review  those 
papers  which  should  be  of  particular  interest  to  one  who 
wishes  to  study  ring  atmospheres. 

General  Theory 

1.  Blamont,  J.   1974.   "The  'Atmosphere'  of  the  Rings  of 
Saturn,"  in  NASA  SP-343,  edited  by  Palluconi  and 
Pettengill,  125-129. 

Blamont 's  paper  is  the  transcript  of  his  comments, 
delivered  at  a  planning  conference  for  the  Pioneer  II 
mission.   His  theoretical  model  is  not  detailed  but 
results  for  the  H2O,  H,  and  OH  atmospheres  are  given. 
Estimates  of  Lyman-a  emissions  are  also  presented. 

2.  Dennefeld,  M.   1974.   "Theoretical  Studies  of  an  Atmos- 
phere around  Saturn's  Rings,"  in  Exploration  of  the 
Planetary  System,  edited  by  Woszczvk  and  Iwaniszewska , 
471-481. 

Dennefeld  was  Blamont 's  student.   His  paper  predicts 
the  density  of  the  ring  atmosphere  (and  that  associated 
with  Titan).   The  sources  considered  are  sublimation, 
bombardment  by  meteoroids,  solar  and  interstellar  wind 
effects,  and  photodisssociation  of  H2O  on  the  ring 
particles.   The  resulting  estimates  are  outdated  for  the 
most  part,  primarily  because  the  production  mechanisms 
are  now  better  understood .   The  conceptual  model  for 
calculation  of  atmospheric  densities,  however,  is  of 
enduring  significance.   It  concisely  gives  a  qualitative 
way  to  describe  the  atmosphere.   One  questionable  asser- 
tion made  in  the  paper  is  that  the  H  atmosphere  will 
be  spherical  due  to  collisions  between  the  H  atoms  and 
ring  particles;  this  is  at  odds  with  the  results  of 
the  present  thesis.   Also,  Dennefeld  estimates  the 
collision  probability  to  be  equal  to  the  optical  . 
thickness . 
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Ballistic  Transport  Models 

1.  Durisen,  R.H.   1984.   "Transport  Effects  due  to  Particle 
Erosion  Mechanisms,"  in  Planetary  Rings,  edited  by 
Greenberg  and  Brahic,  416-446.   University  of  Arizona 
Press . 

Durisen  presents  a  detailed  model  aimed  primarily  at 
prediction  of  "chip  halo,"  i.e.,  macroscopic  ejecta 
resulting  from  meteroid  impact,  formation  (the  model  is 
in  general  applicable  to  the  ring  atmosphere,  too). 
The  assumption  of  "prompt  absorption"  (i.e.,  particle- 
molecule  collisions  are  inelastic)  is  made  for  calcu- 
lations in  his  paper.   In  calculating  collision 
probability,  Durisen  accounts  for  "slant  path"  of  a 
molecule  through  the  rings,  but  it  is  not  clear  whether 
he  considered  the  relative  speed  of  the  particles  versus 
the  molecule.   The  results  of  his  studies  indicate  that 
transport  of  eroded  material  across  the  rings  may  be 
significant  for  understanding  the  fine  structure 
observed  in  the  rings. 

2.  Ip,  W-H.   1983.   "Collisional  Interaction  of  Ring 
Particles:   the  Ballistic  Transport  Process,"  Icarus 
54:  253-262. 

Ip  develops  a  transport  model  independently  of  Durisen. 
The  results  of  the  two  different  approaches  are  quite 
compatible.   Ip  is  able  to  assume  partially  elastic 
collisions,  since  he  employs  the  Monte  Carlo  method  in 
his  simulation.   He  takes  the  collision  probability  to 
depend  on  optical  thickness  only;  a  more  detailed 
estimate  of  Pc  could  easily  be  incorporated  into  his 
model,  though. 


Sources  (Production  Mechanisms) 


1.   Carlson,  R.W.   1980.   "Photosputtering  of  Ice  and  Hydrogen 
around  Saturn's  Rings,"  Nature  2  83:  4  61. 

Carlson  argues  that  photodissociation  of  H2O  on  ring 
particles  could  account  for  the  relatively  high  popu- 
lation of  H  atoms  actually  measured  in  the  vicinity 
of  the  A  and  B  rings.   Note  that  he  calculates  collision 
probability  using  l-e~T  with  x  equal  to  one.   He  takes 
the  sticking  coefficient  to  be  0.22,  but  this  value  is 
based  on  laboratory  experiments  (Brackmann  and  Fite,  1961) 
which  do  not  closely  match  the  conditions  in  Saturn's 
rings.   The  overall  influence  of  the  uncertainty  in 
these  two  factors  should  not  affect  the  order  of  magni- 
tude of  his  calculations,  but  probably  limits  the 
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precision  of  his  estimation  of  total  production  rate 
needed  to  produce  the  observed  H  atom  population. 

2.  Cheng,  A.F.  and  Lanzerotti,  L.J.   1978.   "Ice  Sputtering 
by  Radiation  Belt  Protons  and  the  Rings  of  Saturn  and 
Uranus,"  Journal  of  Geophysical  Research  83:  2597-2602. 

This  paper  is  a  detailed  summary  of  sputtering  effects 
on  ring  particles. 

3.  Cheng,  A.F.,  Lanzerotti,  L.J.,  and  Pironelle,  V.   1982. 
"Charged  Particle  Sputtering  of  Ice  Surfaces  in  Saturn's 
Magnetosphere,"  Journal  of  Geophysical  Research  87: 
4567-4570. 

Using  information  derived  from  the  Pioneer  II,  Voyager  I 
and  Voyager  II  missions,  the  authors  reconsider  the 
sputtering  mechanism,  revising  the  result  of  their 
earlier  (1978)  paper.   Since  the  actual  ion  fluxes 
observed  were  lower  than  those  they  assumed  in  1978, 
the  authors  conclude  that  ion  sputtering  probably  does 
not  contribute  significantly  to  the  H  atom  atmosphere 
in  the  vicinity  of  the  A  and  B  rings  of  Saturn. 

4.  Ip,  W-H.   1978.   "On  the  Lyman-alpha  Emission  from  the 
Vicinity  of  Saturn's  Rings,"  Astronomy  and  Astrophysics 
70:  435-437. 

This  short  paper  suggests  that  H+  ions  escaping  from 
Saturn's  atmosphere  may  collide  with  ring  particles 
where  they  are  neutralized  and  re-emitted,  thus 
contributing  to  the  neutral  H  atmosphere  in  the  vicinity 
of  the  rings. 

5.  Morfill,  G.E.  et  al.   1983.   "Some  Consequences  of 
Meteroid  Impacts  on  Saturn's  Rings,"  Icarus  55:  439-447. 

Morfill  et  al .  make  careful  estimates  of  the  meteroid 
flux  in  the  vicinity  of  Saturn.   Using  this  estimate, 
coupled  with  a  model  for  analyzing  the  effects  of  meteroid 
impact  on  ring  particles,  they  proceed  to  make  predictions 
about  the  ring  halo,  plamsa  production,  and  the  ring 
atmosphere.   The  authors  conclude  that  meteroid  impacts 
should  contribute  significantly  to  the  ring  atmosphere. 


Ring  Characteristics 

1.   Cuzzi,  J.  et  al.   1984.   "Saturn's  Rings:   Properties  and 
Processes,"  in  Planetary  Rings,  edited  by  Greenberg  and 
Brahic,  75-199.   University  of  Arizona  Press. 

This  paper  is  an  up-to-date  summary  of  information  on 
Saturn's  rings.   It  contains  a  short  section  devoted 
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to  the  ring  atmosphere,  which  concisely  explains 
current  knowledge  on  that  subject.   Also,  a  special 
interest  is  an  appendix  which  outlines  radiative  trans- 
fer theory  (which  is  so  important  in  deducing  ring 
characteristics  from  observations). 

2.  Esposito,  L.W.  et  al .   1984.   "Saturn's  Rings:   Struc- 
ture, Dynamics  and  Particle  Properties,"  in  Saturn , 
edited  by  T.  Gehrels  and  M.  Matthews,  463-545.   Univer- 
sity of  Arizona  Press. 

This  paper  nicely  complements  the  preceding  paper 
(Cuzzi  et  al.  1984) . 

3.  Marouf,  E.A.  et  al.   1983.   "Particle  Size  Distributions 
in  Saturn's  Rings  from  Voyager  I  Radio  Occultation, " 
Icarus  54:  189-211. 

This  detailed  paper  provides  the  first  precise  estimates 
of  particle  size  distributions  ever  given  for  the  rings 
(results  are  calculated  for  four  locations  in  the  rings) . 
The  theory  section  in  the  paper  also  provides  a  clear 
application  of  radiative  trasnfer  theory. 

4.  Weidenschilling,  S.  et  al.   1984.   "Ring  Particles: 
Collisional  Interaction  and  Physical  Nature,"  in  Planetary 
Rings ,  edited  by  Greenberg  and  Brahic ,  367-415.   Univer- 
sity of  Arizona  Press. 

Besides  providing  an  in-depth  analysis  of  the  effect 
of  collisions  between  ring  particles  on  the  evolution 
of  the  rings,  this  paper  gives  information  useful  in 
the  study  of  particle  properties. 


Miscellaneous 

1.  Chandrasekhar ,  S.   1960.   Radiative  Transfer.   Dover 
Publications,  Inc. 

Chandrasekhar  produced  the  definitive  work  in  the 
field  of  radiative  transfer  (radiative  transfer  theory 
seeks  to  specify  the  radiation  field  produced  when 
incident  radiation  is  scattered  and  absorbed  by  media 
consisting  of  discrete  particles).   This  book  is  rigorous 
in  approach,  in  fact  it  can  be  appreciated  as  a  work  in 
mathematical  physics,  alone.   It  is  not  light  reading. 

2.  Pollack,  J.B.   1975.   "The  Rings  of  Saturn,"  Space 
Science  Review  18:  3-93. 

Pollack's  paper  is  a  review  of  knowledge  on  Saturn's 
rings  up  to  197  5.   Much  has  been  discovered  since  then, 
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rendering  invalid  some  of  the  predictions  in  the  paper. 
This  does  not  distract  from  the  value  of  his  work  as 
a  historical  review.   Reading  this  paper  is  an  inspira- 
tional experience,  because  Pollack  clearly  demonstrates 
how  knowledge  about  the  rings  increased  through  appli- 
cation of  the  scientific  method. 

3.   van  de  Hulst,  H.C.   1981.   Light  Scattering  by  Small 
Particles .   Dover  Publications,  Inc. 

This  book  is  very  readable.   It  essentially  connects 
the  mathematical  physics  of  Chandrasekhar  (19  60)  with 
the  world  of  applied  research.   It  provides  the  means 
by  which  the  measurement  of  light  scattering  can  be 
used  to  deduce  the  properties  of  the  scattering  medium. 
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